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ABSTRACT 

We analyze the X-ray variabiUty of 264 Sloan Digital Sky Survey spectro- 
scopic quasars using the Chandra public archive. This data set consists of quasars 
with spectroscopic redshifts out to 2; 5 and covers rest-frame time scales up 
to Atgys ~2000 d, with 3 or more X-ray observations available for 82 quasars. 
It therefore samples longer time scales and higher luminosities than previous 
large-scale analyses of AGN variablity. We find significant (^3cr) variation in 
~30% of the quasars overall; the fraction of sources with detected variability 
increases strongly with the number of available source counts up to ~70% for 
sources with >1000 counts per epoch. Assuming the distribution of fractional 
variation is Gaussian, its standard deviation is ~16% on >1 week time scales, 
which is not enough to explain the observed scatter in quasar X-ray-to-optical 
flux ratios as due to variability alone. We find no evidence in our sample that 
quasars are more variable at higher redshifts {z > 2), as has been suggested 
in previous studies. Quasar X-ray spectra vary similarly to some local Seyfert 
AGN in that they steepen as they brighten, with evidence for a constant, hard 
spectral component that is more prominent in fainter stages. We identify one 
highly-variable Narrow Line Seyfert 1-type spectroscopic quasar in the Chandra 
Deep Field-North. We constrain the rate of kilosecond-timescale flares in the 
quasar population using ^8 months of total exposure and also constrain the dis- 
tribution of variation amplitudes between exposures; extreme changes (> 100%) 
are quite rare, while variation at the 25% level occurs in <25% of observations. 
[O III] A5007A emission may be stronger in sources that are identified as X-ray 
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variable; if confirmed, this would represent an additional link between small-scale 
(corona) and large-scale (narrow line region) AGN properties. 

Subject headings: galaxies: active — galaxies: nuclei — X-rays: general 



1. INTRODUCTION 

The variation observed in spectra of active galactic nuclei (AGN) is governed by physical 
processes that we do not fully understand, but are integral to disk/corona emission and 
absorbing outflows. The short time scales over which variation can be observed indicate 
that, in many cases, it is occurring in relatively small structures near to the supermassive 
black hole (SMBH). AGN variability studies provide new temporal constraints for accretion 
and outflow models, particularly concerning size scales; black hole masses and accretion rates; 
ionization structure; emission mechanisms; and relations between various structures such as 
the disk, corona, jet, broad line region (BLR), and narrow line region (NLR). Variability 
studies also attempt to map out how accretion and outflows may depend on luminosity, 
black hole mass, accretion rate, and redshift, in order to identify the fundamental factors 
that influence AGN structure, black hole growth, and galaxy evolution. 

The AGN i dentified in optical surveys such as the Sloan Digital Sky Survey (SDSS; 



York et al.ll2000[ ) can be efficiently studied in other wavebands using publicly-archived data 
from observatories such as Chandra. Because the SDSS spectroscopic quasar catalog is 
primarily optically-selected, the X-ray properties of these AGN have had little influence on 
the AGN selection process. High-quality SDSS spectra also provide secure redshifts, allow 
us to distinguish highly-absorbed broad absorption line (BAL) AGN, and even support 
estimates of bl ack hole masses. Fur thermore, the flux limits of the SDSS spectroscopic 



quasar catalog f lSchneider et al.ll2010[ ) are well-matched to arc hived Chandra obse rvations 



with a very high X-ray detection rate for non-BAL AGN (e.g., iGibson et al.ll2008al ). 



In this study, we describe the X-ray variabihty properties of hundreds of SDSS spec- 
troscopic quasars that have been observed multiple times by Chandra. Most of the X-ray 
observations were serendipitous, in the sense that the quasar was not targeted by Chandra. 
As a result, the sample is also relatively free of biases that could arise from selecting sources 
for targeting in the X-rays. Although most sources were observed only two or three times, 
the large sample size permits us to quantify variability in ensembles out to redshift z ^ 'b 
and covering time scales of tens-of-minutes to years. 



One of the primary goals of our analysis is to characterize quasar X-ray variability as 
a function of time scale, redshift, luminosity, and optical spectral properties. We measure 



- 3 - 



the extent of X-ray variation to determine whether the scatter in X-ray-to-optical flux ratios 



Gibson et al. 



2008a 



Vagnetti et al. 



is do minated by variation or is largely intrinsic (e.g 
20101) ■ We also provide a new test of claims that variability increases at higher redshifts 
( lAlmaini et al.ll2000l : [Manners et al.ll2002t iPaolillo et al.ll2004l ). and examine its dependence 
on luminosity and black hole mass. Another goal of this study is to compare the variation 
of high-luminosity quasars to current models derived from intensive monitoring of local 
Seyfert AGN. Using hardness ratios and model flts, we determine how spectral shape evolves 
as sources brighten and fade, and test simple models of spectral variation against the data. 
Finally, we examine optical spectra for features that differ between ensembles of signiflcantly- 
variable and non- variable sources. 

Because the most detailed studies of variation in individual AGN have been conducted 
by monitoring local Seyfert AGN, we will be drawing from that body of research to guide 
our current study and to suggest directions for future analyses. X-ray variation time scales 
in local Seyfert AGN are observed to depend on black hole mass and accretion rate (e.g., 
O'Neill et al.l l2005l : iMcHardyi 120101 ) . The correlation between X-ray spectral steepness and 
brightness has been p reviously modeled in detail (e.g.. lTaylor et al.ll2003l : IVaughan fc Fabian 
2004 : iMcHardyl l2010l ). The relationship between X-ray and optical variability is complex, 
with at least two different mechanisms modulating the emission. On shorter time scales 
(days), optical emission can lag X-ray emission, suggesting that variability is partly caused by 
X-rays being reprocessed down to lower energies. Variation on longer time scales (years) may 



be driven by changes in the accretion rate (e.g.. lArevalo et al.ll2008l . |2009[ ). The complexity 
of the temporal relationship between optical and X-ray emission may be due to differences 
in the geometry of the regions that emit at these wavelengths. For example, an AGN with a 
larger black hole mass (Mbh) will generally have a cooler accretion disk, so that the region 
of the disk that emits in the optical may be relatively closer to, and subtend a greater solid 
angle of, a central X-ray emitting corona. Such mass-dependent geometric effects could, for 
example, influence the effectiveness of reprocessing and Compton scattering (e.g.. iMcHardy 
2010l ). This example demonstrates one scenario in which variability studies could determine 
that high-luminosity quasars are not simply a "scaled-up" version of local Seyfert AGN. 

In contrast to Seyfert AGN, much of our understanding of the temporal X-ray properties 
of distant, luminous quasars is derived from lower-sensitivity sky surveys and a limited 
number of resource-intensive targeted observation s. Previous X-r a y vari ability analyses using 
ROSAT observations of radio-quiet quasars by lAlmaini et al.l (120001 . hereafter AOO) and 
Manners et al.l (120021. hereaf t er M02), as wel l as an analysis of Chand ra Deep Field-South 
fCDFS: lGiacconi et al.llioOll : ILuo et al.l[2008l : IXue et al.lEoili ) AGN bv lPaolillo et all kmA 
hereafter Pao04), have found an intriguing tendency for AGN at higher redshifts {z > 2) to 
have larger X-ray variability amplitudes than would be expected from an extrapolation of 
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the properties of lower-redshi ft AGN. A related stu dy using XMM-Newton observations of 
the Lockman Hole region by fjPapadakis et al.ll2008l . hereafter Pap08) found that variability 
decreased with increasing redshift in the sample overall, but for a given luminosity range 
of AGNs, variability increased out to z ~ 1, then remained constant at higher redshifts. 
Pao04 also estimated that a large fraction (> 90%) of AGNs likely exhibit X-ray variability, 
and sources for which spectral variability could be measured exhibited a tendency to soften 
spectrally as they brighten. 

Our current sample allows us to expand on these previous studies in several ways. It is 
constructed using high-quality optical spectra that can be used to unambiguously identify 
quasars and determine their redshifts. By contrast, AGN were estimated to account for only 
~80% of the CDFS sample of Pao04. Photometric redshifts were used for the CDFS sources, 
and some of the less-luminous {Lx < 10^^ erg s~^) sources may have been contaminated by 
emission from their host galaxies. Drawing on the large area of sky covered by the SDSS 
survey, our sample includes a large number of highly- luminous quasars, and extends to 
luminosity levels ~10 times higher than even the AOO and MOO ROSAT samples. (See 
§2.41 for further description of our sample properties.) Compared to earlier i?0S'j4r studies, 
Chandra^s sensitivity to hard X-rays permits us to measure hardness ratios and spectral 
shapes at energies >2 keV, where absorbing material (if present) has a weaker effect on 
spectral shape. Chandra's spatial resolution also resolves away background contaminants. 
The data in our serendipitous sample cover long time scales, w ith rest-frame time s between 
epochs up to 5.4 yr. By contrast, the Deep ROSAT Survey fIShanks et al.lll99ll ) used for 
AOO and much of the M02 sample spans about two weeks in the observed frame, while the 
Lockman Hole observations of PapOS covered under two months and the CDFS observations 
used by Pao04 were collected over about 15 months. Our sample also has a larger number of 
sources at higher redshift {z > 2), although larger samples are still needed in this regime. For 
these reasons, we especially focus on the dependence of variability on redshift, luminosity, 
and time scale, and also examine how spectral properties are related to X-ray variation. 

Although our current analysis is focused on X-ray variation, w e note that the te mpo- 
ral em phasis of upcomin g deep- wide surveys such as Pan-STARRS flKaiser et al.ll2002[ ) and 
LSST (jivezic et al.ll2008l ) will greatly enhance variability studies by selecting large new sam- 
ples of bright AGN and also by extending our understanding of the temporal properties of 
AGN in optical wavebands. Optical and X-ray views are complementary because the pro 



cesse s that generate optical and X-ray emiss i on are strongly rela ted (e.g., IZamorani et al. 



1981; 


Strateva et al. 


2005; 


Steffen et al. 


2006; 


Gibson et al. 


2008a 



tical data sets will support increasingly sophisticated AGN research that incorporates time- 
domain information across the spectrum. The SDSS has measured two epochs of photometry 
for large numbers of spectroscopic quasars, permitting the construction of ensemble struc- 
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Vanden Berk et al. 


2004) 


Markowitz et al. 


2003h. t] 


(e.g., 


Wilhite et al. 


2005) 



(e.g. 
e.g. 



2003[). t he optical continua of SDSS quasars become bluer as they brighten 



20051 ). As for X-ray studies, the amplitude of fractional variability de- 
creases with increasing luminosity, an d a positive correlation of variability amplitude wit h 
redshift has also been observed (e.g., ICid Fernandes et al.ll2000l : IVanden Berk et al.ll2004l ). 
The optical structure function is well-represented as a po wer law for time s cales of < 1 yr, 
but appears to flatten at longer (multi-year) time scales ( llvezic et al.l 120041 ). Some models 
invoking a combination of random, discrete emi ssion events are disfavo r ed because they do 
not r eproduce optical structure functions (e.g., iKawaguchi et al.l llQQSt I Vanden Berk et al. 



20041 ). Damped random walk (DRW) models fit t o individual AGN light curves generally 
indicate a damping ti me scale of ~100 days or more flKelly et al.ll2009t iKozlowski et al.ll2010l : 



MacLeod et al.ll2010[ ). corresponding to structure function flattening at longer time scales. 



In the following sections, we describe the sample selection and data reduction process 
([|2]), explain the procedures we use to formally detect and characterize X-ray variability ([|3]), 
discuss some physical implications of our results (^, and provide a concluding summary 
(©• Throughout, we use a cosmology in which Hq = 70 km s~^ Mpc~^, Qm = 0.3, and 
= 0.7. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. The Chandra Sample 

In order to obtain a list of quasars observed by Chandra, we searched the Chandra 
archive to find all AC IS-S or ACIS-I observations of SDSS Data Release 5 (DR5) quasars 



( ISchneider et al.l 120071 ) that used no gratings and were public as of 13 Jan 2009. We use 
this sample to design and calibrate our complex data-analysis pipeline for an initial study. 
In a follow-up study, we will expand the sample to include recently-observed quasars and 
additional AGN selected using new metrics such as optical/UV variability. For each spec- 
troscopic quasar, we identified candidate Chandra observations that had telescope aimpoints 
within 15 arcmin of the given quasar. Then, we investigated these candidate observations to 
determine whether a quasar fell on an active detector chip. In cases where a source fell within 
32 pixels of a chip boundary, we discarded the candidate observation because of increased 
uncertainty of the instrument response in those regions. Where necessary, we corrected as- 
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pect offsets using the prescription available onlinejj Our final sample of multiply-observed 
sources includes 763 Chandra observations of 264 SDSS quasars. 

We reduced each observation of an SDSS quasar using CIAO 4.1.20 following the proce- 
dures listed online for the tool psextract Jf] This tool does not handle cases where zero counts 
are present in the source extraction region. We flagged these cases for special handling. For 
each quasar observation, we generated instrument response files (ARFs and RMFs) either 
directly (using the mkarf and mkrmf tools for zero-count sources) or indirectly (through 
psextract). These response files include a correction for the buildup of contaminant on 
the ACIS chips. The instrument response accounts for spatial and temporal variation as a 
function of detector position and time. 

For point sources brighter than r = 20 mag, SDSS astrometry is accurate to ~50 milli- 
arcecondfl or about 10% ot a Ckanira pixef We can therefore use SDSS astrometry to 
localize sources on Chandra CCDs. We performed "forced photometry" on the X-ray images, 
extracting source counts from a circular region with radius equal to the 90% encircled energy 
fraction at 1.5 keV for a given off-axis angle. This was done even in cases where the source 
is not detected in an image. Extraction radii were computed using data tabulated on the 
Chandra X-ray Center web sitejf] 



2.2. Background Estimation 

Backgrounds were generally extracted from an annular region centered on the source 
position. We selected the inner and outer radii of the annulus to be IS + r^ and 45 + rs pixels, 
respectively, where is the source extraction radius. In our experience, this prescription 
provides sufficient area to estimate backgrounds reliably at small off-axis angles without 
becoming large enough to be contaminated by numerous sources at large off-axis angles. 
In cases where background regions fell partly off-chip, we visually selected a new, circular 
region that was near to the source and appeared not to be contaminated by other sources. 
In addition to visual inspection for source contamination, we also checked each background 
programmatically, searching for sources detected by the Chandra tool wavdetect that fell 

ihttp://cxc.harvard.edu/cal/ASPECT/fix_offsct/fix_offset.cgi 

^http: / / cxc.liarvard.edu/ciao/ 

■^http:/ /cxc. harvard.edu/ciao/thrcads/pscxtract/ 

^http://www. sdss.org/dr7/products/general/astronietry.html 

■5http://cxc.harvard.cdu/cal/Hrina/psf/ECF/hrmaD1996-12-20hrcixcLN0002.fits 
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inside or near our background regions. We generated a list of wavdetect sources using a 
conservative sigthresh threshhold of 10~^, which roughly corresponds to 10 false detections 
per chip. In cases where sources were found to contaminate our background regions, we 
selected a nearby background region that was free of contamination. 



2.3. Sample Properties 



Figure [T] shows the distribution of the number of times a source was observed by Chandra 
in our sample. Of 264 sources, 182 were observed 2 times, while the remaining 82 were 
observed 3 or more times. Four sources were observed 15 or more times. Figure [2] shows 
example light curves for these sources. 

Some analyses depend on the time-domain and redshift coverage of the archive data. 
To illustrate this coverage. Figure E] shows the maximal rest-frame time span (max Atgyg) 
between observing epochs for each source in our sample. Time spans between archived 
observations range from hours to >5 yr for lower-redshift sources. For sources with z > 2, 
time spans of up to ~1 yr are covered. Of course, we can also examine time scales shorter 
than the maximum value of Atgys for sources observed more than two times. 



In an earlier study, iGibson et al.l ( l2008bl ) found that Chandra detects nearly all (~100%) 
non-BAL SDSS spectroscopic quasars in observations having exposure times >2.5 ks or off- 
axis angles <10 arcmin. [The Chandra detection rate of SDSS spectroscopic quasars is 
lower, but still very high, for source s with shorter exposures and/or larger off-axis angles up 
to 15 arcmin (IGibson et al.ll2008bl ).] Motivated by these criteria, we define a high-quality 
sample, which we call "Sample HQ", consisting of 167 sources. To construct Sample HQ, 
we culled all observations of any source that had exposures <2.5 ks and off-axis angles 
>10 arcmin. We also required that sample-HQ observations be performed at an off- axis 
angle >1 arcmin, in order to eliminate bias that could be introduced by the target-selection 
criteria of the Chandra observatory. A typical source has about 5 more counts per epoch (64 
vs. 59 counts, on average) in sample HQ than in the full sample. 

For each quasar, we calculate (or estimate) the monochromatic luminosity Lj^(2500 A) 
in units of ergs s~^ Hz~^, which we denote Iv25oo- This value is either calculated directly from 
our fits to the quasar continuu m at (rest-frame) 25 A, o r is extrapolated by normalizing the 
composite quasar spectrum of IVanden Berk et al.1 ( 1200 ll ) to match the observed spectrum in 
the SDSS bandpass . In the former case, continuum fits were performed using the method of 
Gibson et al.l ( 120091 ). in which a reddened power law or a polynomial was fit to each spectrum, 
excluding regions with broad emission lines, BALs, or ionized iron emission. A plot of Iv25oo 
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as a function of redshift is shown in Figure S] for our full sample. The bulk of the sample 
spans a factor of about 40 in luminosity, from lO^^'® < L2500 < 10^^'^ erg s~^ Hz~^. Because 
the SDSS survey is flux-limited, L2500 is strongly correlated with redshift. Throughout this 
work, we use quasar redshifts reported in the SDSS quasar catalog. These software-generated 
redshifts were verified by visual inspec tion. Repeat observ ations of SDSS quasars have shown 
an rms difference in redshift of 0.006 (jWilhite et al.ll2005l ): this level of precision is sufficient 
for our study. 

We also make use of central black hole masses (Mbh) determined from bro ad emission 
line fits to H/3, Mg II, and C IV emission lines for SDSS spectroscopic quasars (IShen et al. 



20081 ). We adopt these values with the understanding that accurate determinati on of Mbh 



i s an a rea of active research, with known discrepancies among existing methods. IShen et al. 



(I2OO8I ) note that measurement errors for the quasar, continua, and emission lines are generally 
dominated by statistical uncertainties (>0.3-0.4 dex) in the calibration of M^j^-estimation 
methods as well as unknown systematic effects in the use of, e.g., C IV emission as an 
estimator. Figure [5] shows estimated black hole masses and redshifts for quasars in our 
sample. For comparison, we have also plotted black hole masses for Seyfert AGN and 



quasars determin e d through reverberat i on mapping and oth er methods (IBian fc Zhadl2003 
Botte et al.1 booi IPeterson et all booi lO'Neill et allboosh . (We note that in some cases 



these masses may be controversial; this plot is simply intended to be illustrative.) Our SDSS 
quasar sample extends well beyond the Seyfert AGN regime to higher black hole masses 
{Mbh > IO^Mq) and, of course, higher redshifts. 

We id entify subsarnples of quasars that are known to be radio-loud, foll owing the 
method of iGibson et al.l J2008bh . Radio fluxes from FIRST jBecker et al.lll995h or NVSS 
(ICondon et al.lll998l ) were used to determine the ratio of flux densities at 5 GHz and 2500 A. 
Sources having a ratio of R* >10 (log(i?*) > 1) were flagged as radio-loud. We classify 
25 quasars as radio-loud by this criterion. For 97 quasars, we are not able to constrain 
log(-R*) < 1 given the approximate 1 mJy (2.5 mJy) limit of the FIRST (NVSS) surveys. 
Almost all (93) of these quasars are at most radio- intermediate (RIQs; 1 < log(-R*) < 2) , 



Miller et al 



2OIII) 



which generally have X-ray properties similar to radio-quiet quasars (e.g., 
at least in single-epoch analyses. Based on the fraction of radio-loud quasars in the subset of 
sources that we can classify unambiguously (i.e., radio sensitivity extends to log(i?*) < 1), 
we roughly anticipate RIQ contaminants in the radio-quiet quasar set. 



We used the SDSS DR5 catalog of Broad Absorption Line quasars (IGibson et al.ll2009l ) 
to identify any sources known to host absorbing BAL outflows along the line of sight. BAL 
quasars typically have strong X-ray absorption, so we treat them separately. We expect some 
contamination from unidentified BAL quasars at lower redshift {z < 1.7) because the C IV 
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BAL-absorption region for these sources does not lie in the SPSS spectral bandpass; BALs 
are identified in ^1 5% of higher redshift AGN fe.g. JHewett fc Foltzll2003l : IXrump et al.lboOG 



Gibson et al.ll2009[ ). In the full sample of 264 sources, we unambiguously identify 18 cases 



of BAL quasars. 



2.4. Count Rates 

For each ACIS chip that observed an SDSS quasar, we construct a weighted exposure 
map using the prescription given at the Chandra X-ray Genter web sitej^ This exposure map 
can be used to estimate the source count fiux from the count rate observed in a given epoch. 
We used a Galactic-absorbed power law for the exposure map weights. In order to roughly 
estimate the shape of the AGIS spectrum obtained for each quasar observation, we fit each 
background-subtracted spectrum with a power law. This mod el was fit t o counts in the 



observed frame 0.5-8 keV energy band using the Gash statistic (IGashl 119791 ). Our goal was 
not to model features in each spectrum, but simply to describe the overall spectral shape in 
terms of a photon index (F) that would be used to construct the weighted exposure maps. 
We constructed exposure maps for photon indices of 1.6, 1.8, 2.0, and 2.2, and selected the 
exposure map corresponding to the value of F that most closely matched our fit value for 
each observation. 

We extracted total counts from a circular aperture and estimated background counts 
from an annular (or offset circular) aperture as described in §2.11 and §2.21 Both total and 
background counts were multiplied by a factor of 1/c^ as an "aperture correction." Because 
the fraction of encircled counts depends on photon energy, we used a different correction 
factor for our soft (0.5 - 2.0 keV), hard (2.0 - 8.0 keV), and full (0.5 - 8.0 keV) bands. 
We calculate the correction factor for each observation by determining the median energy 
of a count in the appropriate band. The correction factor ca does not differ greatly from 
the aperture correction ca = 0.90 corresponding to the 90% encircled fraction for 1.5 keV 
photons for our extraction radii. For the soft, hard, and full bands, typical values of ca are 
0.92, 0.87, and 0.91, respectively. 

Source counts were estimated by subtracting the background from the total number 
of counts. We estimated the la upper an d lowe r limits on the number of source counts 



by applying Equations 7 and 11 of iGehrelsl (119861 ) to the numbers of total and background 



counts, then propagated the errors to determine errors in source count rates. We applied a 
small correction factor to account for the influence of Galactic absorption in each band; the 



''http:/ /cxc. harvard. cdu/ciao/thrcads/spcctraLwcights/indcx. html 
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correction factor was determined using the Galactic column with our power law model fits to 
each spectrum. We then divided the number of counts by the average value of the weighted 
exposure map in the aperture to obtain count rates (in counts cm~^ s~^) in the soft, hard, 
and full bands for each observation of an SDSS quasar. Properties for individual sources 
are given in Table [H while properties measured for individual epochs are given in Table [2] 
and Table |3]0 The count rates used in this work are measured in the observed frame unless 
otherwise specified. Our statistics are often dimensionless, so if we apply a factor of (1 + z) 
to the time dimension when we calculate count rates or derived quantities such as count rate 
errors, this factor would cancel out of the final statistic. See, e.g., equations |2l HI and [51 For 
many of our analyses, we do not classify sources as "detected" or "not-detected" (at some 
confidence level), but instead work directly with the background-subtracted count rates. For 
faint sources, these rates can be zero or even negative. In any case, there is a high detection 
rate in our sample. For example, full-band detections were obtained at >95% confidence for 
690 of 763 observations in the full sample and 490 of 507 observations meeting Sample HQ 
requirements. Furthermore, for some analyses, we restrict the sample to sources with higher 
numbers of counts in order to maximize signal-to- noise. 

We calculate "count-rate luminosities" from 0.5-8 keV count rates as: 



L, = PAI + zYAtiD^ 



Li 



where Pi is the flux in counts cm~^ s~^ (calculated by dividing net source counts by the 
exposure map), is the luminosity distance, and the exponent x incorporates both a 
bandpass correction and a i^-correction to account for the fact that the Chandra bandpass 
is sampling different spectral regions as a function of redshift. We h ave assumed a power law 



with photon index of F = 1.8 for this correction, so that x = 0.8 (jSchmidt &: Greed Il986l ). 
We use a fiducial photon index of 1.8 instead of our best-fit photon indices because we do 
not want to amplify any effects that are purely due to measurement error or modeling (e.g., 
absorption features that change dramatically between observations). We do not convert 
count-rate luminosities into traditional luminosities (in units of erg s~^) because doing so 
would require assumptions that decrease the accuracy of our results. Calculating the typical 
photon energy used as a conversion factor in a given bandpass involves making an assumption 
about the spectral shape or performing a model fit; both of these can introduce additional 
errors, especially for atypical and/or faint sources. 

Figure [6] compares the luminosity and redshift distribution of our sample to that of AOO, 



^In these tables, we list some measurements, such as count rates and their errors up to three digits after 
the decimal. Of course, the measurements should not be considered significant to this number of digits. This 
approach is adopted to avoid accumulating round-off errors and assist machine interpretation of the data. 
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Pao04, and Pap08. (Data for this figure were kindly provided by O. Almaini and M. Paolillo; 
the data for the M02 sample were not available.) For illustrative purposes, we estimate 
Lo.5_8, the luminosity in the 0.5-8.0 keV band, from the count-rate luminosities in our sample 
assuming that the energy of each photon is ~1 keV. The Pao04 sample drawn from the CDFS 
and the PapOS Lockman Hole sample cover lower luminosities (and shorter time scales), while 
our sample approximately covers the AOO sample and extends it to luminosities up to ~10 
times higher. Our sample includes more quasars at higher redshifts (with 63 quasars at 
2 > 2), as well. Figure [7] shows the maximum rest-frame time scales probed by each study, 
estimated for sources in AOO, Pao04, PapOS by dividing the maximum observed-frame time 
scale by {1 + z). Using the Chandra archive, we can probe time scales up to 10 times longer 
than in the lower- luminosity sample of Pao04, and ~100 times longer than in AOO and PapOS. 
Figure [7] shows (as a gray box) power spectrum break time scales estimated for some typical 
sources in our sample using the relation given in lMcHardy et al.l ( 120061 ). The X-ray archives 
enable us to explore the regime beyond this break time scale. 



3. ANALYSIS 

3.1. Identifying Variation 

For each quasar, we performed a one-parameter fit to the photon fluxes using our esti- 
mated errors to determine the best constant flux that would describe our data. As expected, 
we find that the observed distribution of xl is skewed to higher values than would be ex- 
pected if the constant model were a good fit to the data. Of course, several factors affect 
the distribution of xl in our analysis, including non-normality of the count rate distribution 
and the estimation of source count rate errors. 

For reasons such as these, we adopt a different method to determine whether sources 
can be robustly classified as "variable." For each epoch of a given source, we calculate 
the expected total (source -|- background) count rate in our extraction aperture from the 
measured background and best-fit constant fiux. We fiag an epoch as "variable" if the 
observed count rate is higher or lower than the number of counts corresponding to a deviation 
from the best-fit value at >99% confidence, according to a Poisson statistic. Any source with 
at least one variable epoch is considered a "variable source." For the full sample, 74 of 264 
sources are classified as variable; 54 of 167 sources are variable in Sample HQ. We classify the 
remaining sources as "non-variable," with the caution that this term should be understood 
to mean that we did not detect variability within the limits of our data. 
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3.2. Sensitivity to Variation 

Variability tests are more sensitive for sources witli larger numbers of counts. In order to 
examine tliis dependency in our data, we plot in Figure |8]the mean number of source counts as 
a function of redshift. The "mean counts" are calculated for each source using only epochs 
from Sample HQ that have high exposures and low off-axis angles. Of course, the mean 
number of counts depends on variability characteristics, so this figure should be considered 
to be only a rough exploration of our data set. Sources flagged as variable are plotted in 
red, while non-variable sources are plotted with black points. Sources known to host BALs 
are plotted as circles, while known radio-loud sources are triangles and known radio-loud 
BAL quasars are stars. The remaining sources, including those for which radio-intermediate 
status and BAL absorption could not be ruled out ( §2.3p . are plotted as squares. As a rough 
guide, we have also plotted solid curves indicating typical numbers of counts for quasars 
with 0.5-8 keV luminosities of 10^^, 10^^, and 10^^ erg s~^. The curves were constructed 
assuming a hypothetical source with an unabsorbed, F = 2 power-law spectrum observed 
on-axis on an ACIS-I chip for 18 ks (which is a typical, median exposure time in our sample). 
Individual sources may differ from this hypothetical source in various respects. 

Variability is primarily detected in quasars having >50 counts per observation (on av- 
erage), with about half of those sources flagged as variable. Our sensitivity to variability is 
diminished at higher redshifts (2; > 2), where the majority of sources have <50 counts. On 
the other hand, we detect variability in 9 of 14 sources with >700 average counts and 7 of 
10 sources with >1000 average counts at redshifts z < 1.2. 

3.3. Fractional Variation 

In Figure [HI we show the maximum variation from the best-fit constant count rate for 
the sources in Sample HQ as a function of the mean number of counts per epoch. The y-axis 
for the plotted points is the maximum of |r/ro — 1| for all HQ observations of a source, where 
r is the measured flux and tq is the best-fit constant flux for that source. The solid black 
line shows the "Sa limit" relation [y = 3y/x/x for a mean number of counts x), which is a 
rough approximation of our variability-selection criterion. The thick, red line indicates the 
relation y = f{x), where f{x) is the fraction of sources that are identified as variable in the 
set of sources having mean counts < x. For the entire sample including sources with small 
numbers of counts, ~30% of sources are classified as variable at high significance. 

In the following discussion, we calculate fractional variation using all available measure- 
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ments for each source. We define the fractional variation, F, using the formula: 

F = (c,-c,)/(q + c,), (2) 

where q is the count flux in the later epoch and cj is the flux in the earlier epoch. Mathe- 
matically, F represents half of the full distance between two measurements, (cj — Cj)/2, as 
a fraction of the average value of those two measurements, (q + Cj)/2. We have chosen this 
functional form to be symmetric (up to a sign) in Cj and cj and to roughly represent the 
fractional deviation from some "average" flux. Each measurement of F between two epochs 
is associated with a rest-frame time between measurements. At sys ■ 



3.3.1. Fractional Variation Over Time 

Given the complex nature of measurement errors in our data set, we choose to assume 
that the intrinsic distribution of F is Gaussian. We use that assumption to constrain the 
distribution of fractional variation given the observed values of F and errors in F. The 
standard deviation of the Gaussi an distribution, cr(F), is calculated using the likelihood 



method of iMaccacaro et al.l (Il988l . hereafter M88) in bins of 100 pairs of epochs. The errors 
on (j{F) that are shown in the plot were calculated in the same way. The epoch pairs were 
constructed using all pairs of measurements for each quasar. A quasar that was observed 
times would therefore contribute A^(A^ — l)/2 pairs. 

Figure [10] shows our estimated values of cr(F) as a function of At sys- Each point in 
the figure is placed at an x-coordinate corresponding to the median Atgyg in the bin for 
which (t{F) was calculated. The dot-dashed line represents a linear fit of <t{F) at time scales 
> 5 X 10^ s and is parameterized by: 

a{F) = (0.000 ± 0.012) log(At,j,,) + (0.156 ±0.083). (3) 

At Atsys ^ 1 day, the fractional variation is about 15.6%. At short time scales {Atgys ^ 
5 X 10^ s), there is no significant variation detected above our measurement errors. If much 
more data were available to constrain (t{F) as a function of time, we could map out the gap 
in the current plot where the fractional variation "jumps" from insignificant (on the shortest 
time scales) and breaks to a flatter shape at longer time scales. 

We note that one source, J142052.43 + 525622.4, has so many observations that it 
dominates one of the long time-scale bins. It varies somewhat less than the typical source, 
and if it is removed, the fit model shows a mild, but insignificant increase with longer 
Atsys- Interestingly, J142052.43 + 525622.4 has weaker [O III] emission, similar to the trend 
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observed for non- variable sources ( §3.10p . although it is technically classified as "variable" 
in our sample. 

The red dotted lines in Figure [10] indicate the la and 3a upper limits on a{F) for 
15 epochs of radio-loud, non-BAL quasars with Atgyg > 5 x 10^ s. Because the 3a upper 
limit for radio-loud quasars is only a little larger than the standard deviation of fractional 
variation for radio-quiet quasars, radio-loud quasars appear to be less variable than radio- 
quiet quasars in our sample. If real, this effect could be due to additional, relatively-constant 
jet emission that dilutes the variable X-ray spectrum of the disk corona. 

One shortcoming of this "ensemble approach" to measuring fractional variation is that 
sources with more observations have more infiuence on a{F), because they contribute more 
epochs to the ensemble average. We can also measure quasar variation using a single epoch 
from each source and comparing it to the best-fit photon fiux for that quasar. This approach 
has the advantage of placing each quasar on an equal footing, so that a small number of 
quasars do not dominate the results. However, we do not have enough epochs to map 
out variation in detail as a function of Atgyg, and of course there is no particular "time" 
associated with the best-fit constant fiux. If we simply compare one epoch per source to 
the best fit count rate for that quasar, we find a time-independent standard deviation of 
fractional variation (cj/co — 1, where cq is the model fiux) of about 16.7% ± 2%. This result 
is consistent with values of a{F) found above for time scales >1 week, which is reasonable 
given that such time scales represent the large majority of our data set. 



3.3.2. Symmetry of Variation 
Asyr nmetry in optical light c urves can distinguish among different models of AGN emis- 



sion (e.g. 



Kawaguchi et al.lll998[ ). Here we test whether any asymmetry is evident in the 
X-ray light curves of quasars in Sample HQ. We considered the fractional variation, F, be- 
tween the earliest and latest HQ epochs for each source. We find 87 cases with F > and 
80 cases with F < 0, indicating no strong tendency for quasars to get brighter or fainter over 
the time scale of our sample. 

We also tested subsamples on longer and shorter time scales. For times between epochs 
Atsys > 10'^ s, the split is 51 with F > and 48 with F < 0. For Atsys < 10'^ s, the split is 
36 and 32, respectively. Of course, there may be asymmetries at a level below what we can 
detect with the current data; the measurement errors do add some scatter to the distribution 
of fractional variation. With a larger sample, it would also be possible to test for asymmetry 
on a smaller range of timescales. 
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3.3.3. Extremes of Variation 

While our assumption of a Gaussian distribution of fractional variation permits us to 
characterize the typical extent of variation, it does not describe the behavior of outliers or 
extreme variation events. Figure [TT] shows the fractional variation, F, between pairs of epochs 
in Sample HQ. For this plot, we have omitted any cases where the fractional variation is < la 
from zero, in order to clearly visualize any extreme values of F. Black points show typical 
quasars, while red or green points signify known BAL or radio-loud quasars, respectively. 
As the plot shows, F is nearly always < 50%, and F > 100% is apparently quite rare. 

Extreme variation has been observed in rare cases such as the narrow-line type 1 quasar 
PHL 1092, which de creased in flux by a factor of ~200 over ^ 3.2 yr in the rest frame 



( iMiniutti et al.ll2009[ ). We can use our data set to place limits on the frequency of such 
events. Figure [12] shows upper limits on the rate at which new observations sharing our 
Sample HQ properties would be expected to show a given magnitude of fractional variation. 



The limits are constructed by assembling all values of fractional variation, F, over a given 

sys <5X 10^ S, 5 X 10^ < Atsys 



time frame. The plot considers three time frames, Atgys < 5 x 10^ s, 5 x 10^ < Atgys < 10^ s 



and Atsys > 10 s. The upper limits for these time frames are plotted as a function of 
|F| with black, red, and green curves, respectively. For a given value of |Fo|, we calculated 
the upper limit on the y— axis by determining the maximum intrinsic rate of occurrences 
of values |-F| > |Fo| given the number of observed cases with fractional variation < |Fo|, 
according to a binomial statistic and using 95% confidence limits|§ In cases where a source 
was not detected at >95% confidence in an epoch, we conservatively forced the variation to 
be 100%. This was done for 19 of 1157 epoch pairs. 

We caution that these upper limits depend on subsample sizes and include unmodeled 
variation due to measurement errors. They are therefore indended only to constrain the 
rates at which variation greater than a given rate F occurs. (We could limit the impact of 
measurement error by dropping sources that have low average count rates, but that could 
introduce new biases.) For example, fractional variation F >100% should occur in < 10~^'^ ~ 
4% of observations. Fractional variation F > 25% should occur in < 10~°'^ ~ 25% of 
observations. 



^Thc upper limits we calculate represent limits on the levels of variability that an observer might expect 
to measure, assuming that our data sample is representative of their observations. This empirical approach 
may also include a contribution from exceptional outliers such as flares or flux drops that would not be 
modeled in an ensemble power spectrum. 
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3.4. Excess Variance 



We calculate the excess variance of measured count rates (e.g., iNandra et al.l Il997 



Turner et al.lll999l ). as: 



(^EV 



\,j:[{n.-^^f-ai). (4) 



where rti is the count rate in epoch z, /i is the mean of the values n^, and cTj is the average 
statistical error in the measurement of rij. In order to avoid any bias introduced by having 
different numbers of epochs per source, we calculate aEv using only the earliest and latest 
observation of each quasar. is therefore always 2 in this calculation. The excess variance 
statistic is not ideal for our data set, which consists of sources observed a small number of 
times over diverse time intervals. However, we consider it for comparison to other published 
analyses. 

We limit our analysis to 131 radio-quiet, non-BAL sources that have >50 counts per 
epoch, on average, in order to screen out a large number of sources with small or negative 



^"ev estimate the error on aEv using the formula given in [Turner et al.l (119991 ). Figure 
shows aEv as a function of count-rate luminosity. Sources at z > 2 are plotted in red. At 
lower redshifts (z < 2), 49 of 113 sources have a%y < 0.001; these are plotted as black open 
circles at y = 0.001. Similarly, 10 of 18 higher-redshift {z > 2) AGN have a%y < 0.001 and 
are plotted with red open circles along the bottom of the plot. No clear patterns are visible 
in the plots except a tendency for higher-redshift sources (plotted in red) to have higher 
luminosities, due to the flux-limited nature of the SDSS quasar survey. 

In order to test for additional structure in the data, we calculated the mean of the excess 
variance values in bins of 21 sources each. To calculate this mean, we used only lower-redshift 
sources at z < 2. This mean a^v is plotted in Figure [13] as a thick green line. Error bars 
on the data points (placed at the median count-rate luminosity for each bin) represent the 
estimated error on the mean. The green line shows a general trend of decreasing aj^y with 
luminosity, as is commonly observed (e.g., AOO, M02, Pao04, Pap08, and references therein). 
A filled green circle represents the mean a'^y for 18 sources at higher redshifts z > 2; its 
value is less than 0.001. Dashed red lines in Figure [T3] indicate 2a and 3a upper limits on 
the excess variance for quasars at 2; > 2. While we cannot completely rule out the possibility 
of an increase at high redshifts, the upper limits indicate that any increase would be at most 
very small. 

The extent to which a quasar is measured to vary can depend on the time scales over 
which it is observed. Local Seyfert AGN show increasing variations with time up to at 
least a break time scale. As we discuss in detail in ^4.21 it would be difficult to account 
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for this effect. We do not know the shape of the power spectrum for quasars, and many of 
our observations are hkely beyond the break time scale (Figure [7]). Instead, we explore the 
possibility that the increase in variability could be timescale-dependent by calculating aEv 
using only observations that are separated by a certain range of time scales. For example, 
we constructed a subsample of observations with rest-frame time scales <1 yr. We tried a 
variety of system-frame and observed-frame time scales, and in no case did we find evidence 
for increased variation aX z > 2. 

Using similar methods, we find no significant correlation between a\y and the average 
hardness ratio (defined in §3.6p of a source. This agrees with the result of Pap08, who found 
no correlation between variability amplitude and spectral slope F in their lig ht curves. Pap08 



note that this result disagrees with earlier results (e.g., iGreen et al.l Il993l ) indicating that 
nearby AGN with steeper spectra showed larger amplitude variations. While this may be 
an indication that luminous quasars vary differently than local AGN, we caution that many 
factors affect these results, including the time scales probed, the baseline of spectral shapes 
spanned by a sample, and sensitivity to variation in fainter sources. In following sections 
( §3.61 and §3.71) . we examine the related issue of how spectral shape changes for a single 
source as it varies. 



3.5. Variability Dependence on Physical Properties 

The Chandra archive provides repeat observations of quasars having a wide range of 
luminosities, redshifts, and black hole masses. While individual sources may not have suffi- 
cient observations to permit sensitive tests of variability properties, we can characterize the 
physical dependence of variability more effectively in an ensemble of observations. We place 
each observation of a radio-quiet, non-BAL quasar in Sample HQ on an equal footing in this 
ensemble by considering the quantity 

Sij = {fij — Cj)/aij, (5) 

where fij is the 0.5 — 8 keV count rate observed in epoch i for source j, Cj is the best-fit 
constant count rate for source j over all epochs, and CTjj is the error on fij. In the absence 
of variability, we would expect the distribution of Sij values to be roughly Gaussian with an 
rms of 1. The mean value of Sij would be ^ (by construction) in ideal conditions where 
variation is symmetric and not infiuenced by outliers. This is generally the case except at 
the lowest redshift, which is infiuenced by a few (^ 5) outliers that have higher Sij values. 
We use Sij to test the hypothesis that intrinsic variability can be detected in the data at a 
certain redshift, luminosity level, or black hole mass. Of course, determining the amplitude 
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or detailed pattern of any intrinsic variability would require simulations to determine what 
pattern of Sij values could be expected from a given variability model. 

In Figure [TU we plot Sy as a function of redshift. We have binned the set of Sy 
values into 5 bins and calculated the sample standard deviation in each bin. Red data 
points (placed at the median x-value of each bin) indicate the square root of the unbiased 
sample variance, with error bars estimated as Aa = a/ ^/2(N — 1) for a bin with data 
points. (However, we note that scatter in S^j can be driven by non-Gaussian effects such as 
outliers.) The unbiased sample variance reflects the influence of outliers, some of which have 
been clipped out of the plot range for visibility purposes. A better estimate of the overall 
sam ple variability propert ies can be obtained using the median absolute deviation (MAD; 



M„.on„„.t „im which i= much le. ^en^itivc to st.ongly-vanablc outhe..B The 



e.g. 

green line shows this value, (Tmad calculated for each bin, with errors roughly estimated 
using the same formula as for the red data points. 

As redshift increases, the amount of variability that we detect decreases, ctmad for Sij 
decreases to ~1 as redshift increases, indicating that any intrinsic variability is dominated 
by measurement errors in this analysis at 2 > 2. As Figure M shows, the average number 
of counts per epoch for a given quasar also decreases at higher redshift, so that variability 
becomes harder to detect at z > 2. However, we find the same result — Sij decreasing to 
the noise level {Sij = 1) — even when we limit our analysis to Sample HQ sources that have 
10 to 50 average counts per epoch. We do not, therefore, find any indication that variability 
increases at higher redshifts in this data set. 

Of course, these statements simply express the extent to which we can measure vari- 
ability in the data. Even bright AGN at high redshifts are presumably variable at some 
level, and their variability patterns may be complex. In order to demonstrate how a simple 
pattern of variation would appear in the data, we have created a set of simulated data in 
which we fix a constant amount of variation (10%, 20%, or 30%) in each observation of a 
quasar, cmad for the simulated Sij values are plotted as blue lines in Figure [HI These 
decrease with redshift because measurement errors are generally increasing with redshift, 
but the simulated lines do not fall all the way to Sij = 1, because they do possess intrinsic 
variability. The green line corresponding to real quasars declines more steeply with redshift 
than the blue lines, again indicating that we see no significant increase in variability with 
redshift. 



^Thc MAD for a set of values is defined to be MAD(a;) = mediani{\xi — median{x)\) . I.e., it is the median 
of the absolute value of residuals from the sample median. It is less sensitive to outliers than the sample 
variance is, and can be used to estimate the standard deviation a according to a ~ 1.483xMAD. 
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We find similar results when plotting Sij as a function of luminosity or Mbh- ctmad 
decreases to the noise level in the highest Mbh bins. (Tmad decreases with luminosity as 
well, although there is a small (but insignificant) upturn at the highest luminosity values 
(Figure [T5]) . 

Finally, we have attempted to distentangle the degeneracy between redshift and lumi- 
nosity by binning our data on both quantities. Even with only four redshift bins (0-0.5, 
0.5-1, 1-2, > 2), this pushes the limits of available data, leaving generally 10-35 Sij points 
per (L, z) bin. In any case, we still find a general trend for a mad to decrease with both 
luminosity and redshift, and no evidence of an increase at z > 2. 

Throughout this analysis, we have worked with a variability signal that is measured in 
the ohserved-baiae 0.5 — 8 keV band where Chandra is most sensitive. To the extent that 
AGN variation differs from one part of the spectrum to another, our ability to detect variation 
is also redshift-dependent. (See further discussion of this point in §4.21 ) The instrumental 
response varies strongly over the bandpass, further complicating any attempt to separate 
energy- dependent variation from redshift dependence. For the current work, when we state 
that "We do not detect variability" in some instance, that statement should be understood in 
the context of the data that we have available, using the Chandra observed-frame bandpass 
and instrumental response. 



3.6. Variation in Spectral Shape 

Hardness ratios (HRs) — here defined as the ratio HR = {H — S)/{H + S), where H 
and S are the observed-frame count rates in the hard (2.0-8.0 keV) and soft (0.5-2.0 keV) 
bands, respectively — may be used to describe the overall spectral shape for sources that do 
not have sufficient counts for detailed spectral fitting. We have com puted observed-fra me 



hardness ratios for all observations in Sample HQ, using the method of iPark et al.l fl2006[ ) to 
determine la HR errors. We calculated separate weighted exposure maps for the soft and 
hard bands using the spectral slope obtained from the best fit to the full-band spectrum of 
each source. (Narrower bandpasses and smaller count rates do not permit separately fitting 
the soft and hard band spectra.) 

We examine all epochs in Sample HQ (excluding those from BAL and radio-loud quasars) 
to determine how HR varies with luminosity. Here we put all epochs of all sources on an 
equal footing, scaled by their HR and luminosity in the earliest available epoch. For each 
epoch, we calculate the change in HR and luminosity from the first observation of that 
source. These quantities [HR — HR{t = 0) and L/L{t = 0)] are plotted in Figure [TBI for 
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lower-redshift [z < 2, black squares) and higher-redshift {z > 2, open circles) quasars. There 
is a significant anticorrelation (at >99.99% confidence, according to a Spearman rank cor- 
relation test) between the change in HR and the L/L{t = 0) for quasars at z < 2. The 
correlation is not significant at z > 2, although the sample sensitivity is weakened by larger 
measurement errors for fainter sources as well as the fact that softer photons are shifted 
out of the Chandra bandpass at higher redshifts. Defining AHR = HR — HR{t = 0) and 
F = Lph/Lph{t = 0), we obtain the fit 

AHR = (-0.298 ± 0.030) x log(F) + (0.003 ± 0.004) (6) 

for the full sample including all redshifts. 

We constructed toy physical models of the trend between AHR and L/L{t = 0). The 
change in spectral hardness could, of course, be associated with a single power law model 
that changes both shape and normalization, following the relation described in Equation O 
The goal of our toy models is to determine whether we can reproduce Equation [6] with only 
a single variable parameter. 

We modeled several scenarios, including: 1) a double power law with a constant hard 
component (F = 1.3) added to a variable softer component (F = 2), 2) a constant power law 
with variable neutral absorption, and 3) a constant power law with an absorber that varies in 
ionization level. For each toy model, we attempted to determine a set of fit parameters that 
caused the HR and luminosity to vary together in a way that reproduced the trends observed 
in Figure dni The absorption-dominated models (cases 2 and 3) generally produced a much 
steeper trend of HR with luminosity than we observed. In contrast, the double power-law 
model (case 1) reproduced the observed trend reasonably well for a hard, constant component 
having F = 1.3 and a normalization at 1 keV of about 20% of the variable, soft power-law 
component. This model is shown as a red curve in Figure [TBI The impact of a two-component 
model on observations of high-redshift quasars is discussed further in §4.31 



3.7. Spectral Fits for Bright Sources 

For sources with large numbers of counts, we can fit spectra to determine how spectral 
shapes and features vary as a function of source brightness. We have selected 16 radio- 
quiet, non-BAL quasars having at least 500 counts per epoch, on average, and fit them with 
a spectral model consisting of a Galactic-absorbed power law absorbed by one absorption 
edge. The edge is constrained to have an energy threshhold of 0.739 keV in the rest frame, 
corresponding to absorption from the O VII ion. The optical depth of the edge, r, is allowed 
to vary. We have chosen the edge model to roughly characterize a power law spectrum with 



- 21 - 



ionized absorption using a single parameter. In fact, we constrain r > at the la level in 
only 7 of 87 epochs. We obtain similar results using a neutral, rest-frame absorption model, 
so the choice of this model is not strongly biasing our results. While the single-edge model 
certainly does not reproduce the complex spectral features that an absorber may imprint 
on an emission spectrum, it does provide a basic model of variable absorption when spectra 
do not permit careful fitting of absorption models that have more degrees of freedom. The 
results of the fits are given in Table HI 

Figure [T71 shows the distribution of (AFj^, Tij) pairs from our fit models. We ordered the 
Sample HQ observations of each source chronologically, then calculated Ar^j as the difference 
in photon indices between epoch i + 1 and epoch i for source j. Tij is the photon index from 
the first epoch (i.e., epoch i) in the pair. Two outlier points have been clipped in this figure 
for visual clarity. The plot shows a strong anti-correlation between AF and F at the >99.99% 
confidence level, according to a Spearman rank correlation test. The median value of F in 
this sample is Tm ~ 2.16. (The value F = 2.16 is not necessarily representative of the full 
population; for example, the sources chosen for spectral fitting were selected to have larger 
numbers of X-ray counts.) Spectra tend to steepen when they are flatter than this value, 
and tend to flatten when they are steeper. Figure [TS] shows how the best-fit photon index 
changes with the power law normalization. There is a strong correlation (99.97% confidence, 
according to a Spearman rank correlation test) between AFjj and ANorm^j, where ANorm^j 
is the difference in power law normalization values between epochs i + 1 and i. Overall, 
spectra steepen as they brighten (and/or flatten when getting fainter), consistent with what 
we have already observed using hardness ratios for a larger sample of sources (see Figure US] 
and gMD. 



3.8. Hard Flux and Soft Flux 



Some local Se yfert AGN show a linear relationship between hard- and soft-band X-ray 
count rates 



le.g.. 



Taylor et al.l l2003l ). When the linear relationship is extrapolated to a 



zero soft-band count rate, a significant non-zero hard-band flux remains. This suggests a 
two-component model of Seyfert AGN X-ray spectra, in which a constant hard spectral 
component is augmented by a softer, variable power-law. The non-zero offset in the hard 
band comes from the underlying hard spectral component that remains when the variable 
power law component is absent. The shape of the hard spectral component, determined from 
flux-flux plots i n different X-ray b ands, can resemble the expected emission from cold disk 



reflection (e.g.. iTaylor et al.ll2003l ) 



For each radio-quiet, non-BAL source having two or more observations in Sample HQ, 
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we use the IDL FITEXY routine to generate a linear fit to the observed frame hard-band 
(2-8 keV) count rates as a function of soft-band (0.5-2 keV) count rates. That is, we fit for 
a and h in the equation: 

th = ars + b, (7) 

where and are the hard- and soft-band count rates, respectively. The FITEXY routine 
accounts for errors in both the hard and soft band count rates. Our current data do not 
permit a precise test of whether the fiux-fiux relation is best modeled as linear, although 
visual inspection indicates that a linear model works well in most cases with more than 
5 observations. The linear model is not formally rejected by the fit statistic in 12 of 15 
quasars with > 5 observations. Because many sources have only 2 or 3 observations, our 
goal is not to test the linear model, but instead to determine a typical value for yo in the 
ensemble that would correspond to a constant hard-band offset. 

Under the approximation that the y-intercept values are Gaussian distributed, we used 
the method of M88 to determine y^, the mean value of the y-intercept in these fits. The 
mean ?/- intercepts, corresponding to a constant hard-band fiux component, were greater 
than zero at >99% confidence both for lower-redshift Sample HQ sources (51 at 2 < 1) and 
also for higher- redshift sources (116 at z > 1). Then, to estimate the contribution of this 
hard, constant component to the hard X-ray spectrum, we re-ran the fits, but this time 
normalized the hard-band count rates by the mean hard-band count rate for that source. 
The normalized fractional offset has a value of 0.428 ± 0.140 for sources at z < 1 and 
0.753 ± 0.165 for sources at z > 1. The increase in the fractional offset (although with 
a large error bar) suggests that the the shape of a putative constant component would be 
flatter than the shape of the variable component, so that it contributes an increasing fraction 
of flux at higher energies. Local Seyfert AGN show hard-spectrum fractional contributions 
of 20-40% at energies >2 keV, possibly due to Compton reflection from the accretion disk ; 



these fractions increase at higher energies (e.g.. lTaylor et al.ll2003l : IVaughan fc Fabianll2004l ). 
If a similar model holds in the quasar luminosity regime, we might expect the fractional 
contribution of the constant component to increase as higher energies are shifted into the 
Chandra bandpass. (See §4.31 for further discussion about the effects on observations of 
high-redshift quasars.) 



3.9. Intra-Observation Variability 



The counts in Chandra event lists are tagged with arrival times, permitting a search for 
variability in the light curves of each observation. We use a Kolmogorov-Smirnov (KS) test 
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to identify any light curve for which the count arrival times are significantly non-uniform. 
(Here we define "light curve" to mean the time sequence of all counts that fall inside the 
source extraction region.) At a 99% KS-test confidence level, we identify 17 potentially vari- 
able light curves. We have removed five additional observations from consideration because 
visual inspection of their light curves indicates that instrumental factors such as background 
flaring are dominating the light curves. A time-dependent analysis of background variation 
is beyond the scope of this study, but we note that in our bright sample (with total count 
rates > 0.01 s~^, described below), the estimated background count rate is <10% of the total 
rate in ^86% of the light curve exposure. Background effects should not strongly affect our 
overall results. 

Six of the 17 light curves that the KS test flags as variable belong to a single source, 
J123800.91 + 621336.0 (hereafter "J1238") at z = 0.44 with absolute z-magnitude Mi = 
—23. 04. This source was ob served many times as it falls in the Chandra Deep Field-North 



(e.g. 



Alexander et al.l 120031 ). Given that we have searched ^800 light curves for a result at 
99% confidence, we cannot consider the variability criterion to be highly significant for the 
remaining 11 light curves. This result is consistent with our previous observation that the 
level of intrinsic variation between observations falls below our detection limit for time scales 
shorter than ~ 10^ s (^33]). 

Figure [li shows the 4000-5200A region of the SDSS spectrum of J1238 (in black) with 



the quasar composite of IVanden Berk et al.l ( 1200 ll ) overplotted in red for comparison. The 
FWHM of the and Mg II A2800 (not shown) lines is ?»1860 km s-\ while the [O III] 
A5007 line is quite weak in comparison to the composite value. J1238 also shows an excess 
of ionized iron emission in the 4500-4600A region. Fitting a Galactic-absorbed power law 
to the observed-frame 0.5-8 keV region, we find that the spectrum is steeper than typical 
quasars, with 2.2 < F < 2.4 in most epochs (Figure I2U]) . In several epochs, F deviates from 
the median fit value of F = 2.6 by 3.3-5.4(T, including two epochs with flatter spectra (F < 2) 



Pogge 


2000; 


Komossa 


2007) 



AGN are often observed to be highly variable in X-rays, and ar e thought to have smaller 
black hole masses an d correspondingly high accretion rates (e.g., 
19991 : iKomossal 120071 . and references therein) . 



Boiler et al.lll996l : iLeighlv 



J1238 was not detected in the VLA FIRST survey, and is classified as radio-quie t in ou r 
study. It was detected as an unresolved source in the deep radio survey of iRichardsl ( l2000l ). 
with a 1.4 GHz flux of 190 ± 13 fiJj, corresponding to log{R*) ^ -0.004. Although J1238 
does show radio emission at faint levels, it is far below the radio-loud limit for our sample. 



As a second test for short-term variability, we implemented an algorithm to search 
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for flares or dramatic, short-term absorption (or dimming; ) events similar to ev ents that 
have been previously reported at quasar luminosities fe.g.. iRemillard et al.lll99ll ). or from 



the Galactic center at much lower luminosities (e.g.. iBaganoff et al.ll200ll ). The algorithm 
breaks each light curve into segments 1000 s wide. This segment size was chosen to provide a 
reasonable balance between the number of segments per light curve and the number of counts 
per segment. The first stage of our algorithm iteratively determines the baseline count rate 
of the light curve, excluding any candidate flare regions. In each light curve, we identify the 
segment that shows the most significant deviation from the mean count rate (according to a 
Poisson statistic), where the latter is calculated using segments that have not already been 
flagged as potentially variable and are not the current segment under test. If the deviation 
is significant at >99% confidence according to a Poisson statistic, we flag that segment as 
potentially variable. Then we iterate the process to find the next-most- variable segment, and 
so on. When no more segments can be flagged as potentially variable, we have determined 
a "clean" estimate of the baseline count rate of the source. Next, the second stage of our 
algorithm uses this new baseline count rate to determine which segments should be classified 
as variable at >99% confidence. 

Our algorithm identifies hundreds of candidate variable segments out of about 21652 
segments, or 0.69 yr of combined light curves. To evaluate the significance of this result, 
we re-run our algorithm on simulated light curves that have the same mean count rates and 
no intrinsic variability. For each segment flagged as variable in either the real or simulated 
light curves, we record the baseline (expected) count rate and the count rate "fractional 
deviation" F^, defined as: 

= iCots-{C))/{C), (8) 

where Cobs is the observed count rate in that segment and (C) is the mean count rate for 
the light curve. 

The distribution of Fc depends on (C), because weaker variability can be detected at 
higher count rates. In Figure [HI we show the distribution of F^ for count rates higher than 
0.01 counts s~^. We also consider the case of lower count rates separately (not shown). For 
a baseline count rate of 0.01 counts s~^, a significant "event" would have a flux increase of at 
least about 80%, or a decrease of about 70%. For higher baseline count rates, the amplitude 
would be smaller. Events in t he high-count-ra t e sub sample could therefore be relatively 
similar to the flare reported by IRemillard et al.l ( 1l99ll ) in the quasar PKS 0558-504, which 
increased by up to 67% over 3 minutes and laste d for 10-20 minutes. Additional rapid flu x 



changes have been observed for this source (e.g., iWang et al.ll2001t iBrinkmann et al.ll2004l ). 
A more extreme example would be PH L 1092, which was o bserved to brighten by a factor 
of ~3.8 over a rest-frame time <3.6 ks flBrandt et al.lll999[ ). Rapid, strong flares have also 



been observed in P DS 456 flReeves et al 1120001): NLSl-tvpe AGN such as RX J1 702.5+3247 
dciiozzi et allboOlh . IRAS 13224-3809 feoller et al1ll997t lOewangan et al.ll2002[ ): and other 
AGN. 



A KS test finds no evidence in most cases of a significant difference in the distribution of 
Fc between the real and simulated data sets either for the distribution as a whole or for the 
separate cases of absorption/dimming (F^ < 0) and flaring [F^ > 0) events. The difference 
in distributions is signiflcant (at >99.9% confldence) for the absorption/dimming {Fc < 0) 
events in the low-counts case; presumably this is because our simulations have not modeled 
the background and therefore do not reach a "floor" where the counts in a time segment 
are dominated by the background. Based on these results from the KS test, we conclude 
that the distributions of F^ do not differ remarkably from distributions in the null case of no 
intrinsic flaring or dimming. 

Of course, it may still be the case that we are detecting variable segments of light curves 
at a higher rate than predicted from the simulations, even if the distribution of Fc is not 
signiflcantly different. Of the four cases under test (low-count-rate absorption/dimming, 
low-count-rate flaring, high-count-rate absorption/dimming, high-count-rate flaring), two 
show event rates that are marginally higher than expected from the simulations. According 
to a binomial statistic, the probability of seeing more absorption/dimming events in the 
high-counts case is 1.4%, while the probability of observing more flaring events in the low- 
counts case is 3.5%. In each of the four cases, the flare rate is consistent with zero at 99% 
confldence. Upper limits on the rates of 1 ks events (after subtracting off the simulated rates, 
in numbers per observed-frame year) are < 6.7 yr~^ (low-count-rate absorption/dimming), 

< 29.3 yr~^ (low-count-rate flaring), < 40.5 yr~^ (high-count-rate absorption/dimming), 

< 9.0 yr~^ (high-count-rate flaring). The limits quoted are for 99% (single-sided) confldence 
limits. The upper limits are higher in cases that showed marginal evidence for signiflcant 
event rates above the simulated values because a higher (but not highly-signiflcant) number of 
events was detected in these cases. For the full sample (combining both low and high count 
rates), we have the following rates: < 37.3 yr~^ (absorption/dimming) and < 51.5 yr~^ 
(flaring). At most a small fraction of time (< 0.2%) is spent by quasars in the short (ks-time 
scale) flaring or absorption/dimming states that our algorithm is designed to detect. 

The algorithm presented here is only a flrst empirical attempt to characterize any pos- 
sible flaring activity in the large sample of archived light curves. One possible enhancement 
would be to search for flares that cover longer time scales. As a preliminary test, we applied 
our algorithm to search for flares in light curve segments 5 ks long (rather than 1 ks). (At 
0.01 counts s~^, a signiflcant event would require an increase by 36% or a decrease by 32%.) 
Our results are qualitatively similar, but the upper limits we can place are much weaker due 
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to the smaller numbers of light curve segments. 



3.10. Optical Spectral Properties 

The SDSS spectra of our sources have a wide range of spectral shapes and signal-to- 
noise values. In order to conduct a basic comparison of the optical spectral properties of 
variable and non- variable sources in the spectral region that contains the H/3 and [O III]A5007 
emission lines, we construct composite spectra of all sources with redshifts < 2; < 0.8 and 
>50 source counts per epoch, on average. We have selected this region of spectrum because 
the H/3 and [O III] emission lines are commonly used to estimate black hole masses and 
accretion rates, and to identify subclasses of AGN such as Narrow Line Seyfert Is. The 
minimum count rate requirement ensures that we can detect variability at the >50% level 
at ~ 3(T confidence ( §3.2p . 

The composite spectra are constructed by calculating the geometric mean of the spectra 
in a given (rest-frame) wavelength bin. The input spectra are normalized to 1 at 5100 A in 
each case; note that the output spectrum depends somewhat on the choice of normalization 
wavelength. We chose to normalize at 5100 A because it is a commonly-used spectral region 
that is less affected by complex emission or absorption features. We omitted from the 
calculation any radio-loud quasars or known hosts of BAL outflows. The variable composite 
was derived from SDSS spectra of 22 sources, while the non-variable composite represents 
15 sources. 

In Figure [221 "we plot the composite spectra in black for variable (top panel) and non- 
variable (bottom panel) sources. As before, "variable" sources are those for which the X-ray 
count rate is inconsistent with a constant value for at least one epoc h. For comparison, 



we ha ve plotted a renormalized SDSS quasar composite spectrum from IVanden Berk et al. 



( 120011 ) in red in each panel. 



We find at most mild differences in the composite spectra for the H/3 A4861 line or 
the broad pattern of ionized Fe emission in the plotted region. However, [O III] emission 
for non- variable sources is noticeably stronger than typical in three lines (4364, 4960, and 
5007 A) where [O III] is clearly present in emission. Visual inspection of the individual 
spectra supports this observation, although there is much diversity. 

To investigate this trend further, we estimated the monochromatic luminosity vLy at 
5100 A and also the equivalent width (EW) of the [O III] A5007 emission line for each 
spectrum. We used a simple linear continuum fit across the H/? and [O III] emission line 
region to estimate the EW. No attempt was made to deblend ionized Fe emission; although we 
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note that the composite spectra do not indicate a strong difference in Fe emission between 
the two subsamples, in general. Figure |23] shows the results of these calculations. While 
both subsamples (variable and non-variable) cover a similar range of luminosities and EWs, 
the non- variable sources are more likely to occupy the high-EW portion of the plot, with 
\og{EW) > 1.3. A KS test finds no evidence of a difference in the luminosity distributions, 
but a suggestive probability (97% confidence) that the EW distributions differ. 

Because we have limited this study to sources with at least 50 counts per epoch on aver- 
age, our classifications of "variable" and "non-variable" are very similar to placing a cut on 
the amplitude of variability. Still, the individual observations do have different sensitivities, 
so the two concepts are not exactly similar. In order to investigate [O III] strength as a 
function of variability amplitude, we define 

S-j = max(|^,,|), (9) 

where Sij is defined in Equation [5] and the maximum is taken over all epochs i for source j. 
Sj represents the maximum amplitude of variation (in units of a) for each of the 37 sources 
we are investigating. In fact, we find that sources with larger values of Sj are generally the 
sources we flagged as variable; the two criteria are nearly equivalent. If we cut the sample 
at the median value S'j = 2.8, then a KS test indicates that the [O III] distributions of 
the two subsamples are incompatible at 96% confidence. If we place the cut closer to the 
division between the "variable" and "non-variable" sources, at S'j = 1.5 to 1.9, then the KS 
test indicates a sample difference at >99.5% confidence. (We find similar results, though 
less significant by a few percent, using a definition in which variability amplitude is defined 
as = maxj(|(/jj — Cj)/cj\), with the best-fit count rate in the denominator instead of the 
measurement error.) As before, we find that there are some indications of a real effect in the 
data, but larger samples are needed to test this result more sensitively. 

It is also possible that some other effects in the data are influencing our "variable" 
and "non-variable" subsamples. If we had enough sources, it would be ideal to control for 
factors such as the number of observations of a given source and the time scales over which 
the sources were observed. With our small sample, this is difficult to do. Although our 
variability criterion is not associated with a specific time scale (but is simply compared 
to the best-fit constant rate), we can investigate what happens if we remove ten cases of 
"variable" sources that were observed more than two times, under the extreme assumption 
that they would not have been flagged as variable from their first two observations alone. In 
this case, the split in [O III] EWs becomes less significant (with a 19% probability of greater 
difference in EW distributions), but this is expected due to the fact that the remaining 
sample is small. Effectively, this cut has removed many "variable" sources without changing 
the cluster of non-variable sources having stronger-than-average [O III] emission. 
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[O III] emission and X-ray variation are two effects that are associated with the well- 
Icnown "Eigenvector 1" description of AGN. Our current study has been partially motivated 
by a desire to determine if and how such factors are actually related, and what the physical 
significance of such a relation would be. We discuss this issue further in §4.41 



DISCUSSION 



4.1. Chandra Sensitivity and Quasar Variation 



Chandra^s ability to detect significant X-ray variation depends strongly on the number 
of counts obtained for a source. Figure |9] demonstrates the sensitivity of Chandra to variation 
in archived observations of SDSS spectroscopic quasars. A minimum of ~10 counts per epoch 
are needed before variability can be detected at the 100% level (corresponding to a source 
doubling or halving its count rate compared to the best-fit constant rate). For sources with 
fewer than 20-30 counts per epoch, we see minimal evidence of variability. 

Overall, variability is detected in ;^30% of the sample, including sources with low num- 
bers of counts. For sources with larger numbers of counts, variability can be measured with 
greater precision. Eleven of 19 sources having >500 mean counts are detected as variable 
at the ^13% level, and 7 of 10 sources having >1000 mean counts are detected as variable 
at the ^10% level. Based on this small sample, we expect that X-ray variability at the 
> 10% level would b e observed in most, if not all, quasars with sufficient numbers of counts. 
Paolillo et al.l (j2004f ) likewise detected variablity in over half the AGN in their analysis of 
the Chandra Deep Field-South, and estimated that intrinsic variability could be present in 
a fraction approaching 90% of their sources. 

The level of fractional variation ^ 16%) we observe certainly contributes to the 

scatter in observed X-ray-to-optical flux ratios for non-BAL SDSS quasars. This amount 
of variation is smaller than that req uired to explain the scatter in X-ray-to-optical spec- 



tral slo pes as due to variability alone fIStrateva et al.ll2005t ISteffen et al.ll2006l : iGibson et al. 
2008a| ). at least over time scales of <2 y r. Using simultaneous UV/X-ray observations 



from XMM-Newton, IVagnetti et al.l (120 lOf ) have also found that the observed scatter in 
X-ray-to-optical flux ratios cannot be explained by variability alone. Some physical differ- 
ences in quasars must also contribute to an intrinsic scatter in X-ray-to-optical flux ratios. 



On time scales shorter than At 



sys 



2 X 10^ s, we do not detect quasar variability. 
At longer time scales, the amount of intrinsic variability is roughly constant or perhaps 
slowly increasing. If we assume that 5 x 10^ s represents a "break" in the quasar luminosity 
structure function, an assumed linear relationship purely between the break time scale (Tg) 
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and Mbh would unde restimate Mbh in our sources by an order of magnitude or more (e.g. 



McHardy et al.l 120051). For luminous quasars, an additional term dependent on bolometric 



luminosity (e.g.. lMcHardy et al.ll2006l ) may be important to increase the predicted black hole 



mass for high- Mbh, high- luminosity quasars. 



4.2. Variability at Higher Redshifts 

Although earlier studies have sugges t ed that variability increases at higher redshifts 



(lAlmaini et al.l l2000l : [Manners et al.l l2002t jPaolillo et al.ll2004l ). we have not been able to 



confirm that result in our data ( §3.4^ §3.5p . In fact, any increase in excess variance at 
z > 2 would be above the ~ 2a upper limit on the excess variance we actually measure 
( §3.4p . We also do not detect any increase in variation at high black hole masses. Increased 
variability at higher redshifts was not detected even when we cut the sample to remove 
sources with lower numbers of counts, thus maximizing the sensitivity to any variability. 
Looking forward, dedicated X-ray observations targeting higher-redshift quasars that are 
matched to lower-redshift quasars in other sample properties would be ideal for a definitive 
test of redshift-dependence. 

We stress that our approach has been to look for an empirical signature of variation 
using the data and observational bandpass that are available to us. Some earlier studies 
(AOO, M02) have generated correction factors for their measured variance that model the 
effects of sampling sources on different time scales. These corrections rely on an assumed 
shape of a quasar power spectrum to boost the observed variation of sources observed over 
shorter times, as quasars tend to vary more over longer time scales. However, several factors 
argue against attempting such corrections to our study. It is not clear what the shape of the 
power spectrum should be over the wide range of time scales in our data, especially since 
many of our data points sample time scales beyond the power spectrum break (see Figure [7]). 
We know even less about this region of the power spectrum for quasars than we do about the 
slope of the quasar power spectrum on shorter time scales. Furthermore, some statistics in 
our study (such as Sij) do not have a well-defined time scale associated with them. For the 
case of the aEv statistic, which we carefully constructed to depend on only one time scale 
per source, we did not find evidence for increased variability at higher redshifts even when 
we limited the range of system- frame or observed- frame time scales ( §3.4p . 

Applying a time dilation correction factor of (1 + z)^^^ (AOO) would have only a small 
(40%) effect even at z = 3, and would still not result in a detection of "increased variability" 
at higher redshifts in our data. Similarly, the study of P04 notes that "no correction for time 
dilation has been applied to this data since it would require us to make a somewhat arbitrary 
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assumption about the power density spectrum (PDS) of our sources...." That assumption 
becomes even more arbitrary for the longer time scales in our data. There is, of course, 
some danger that if we mis-characterise the correction for our sample, we will be creating 
a redshift-dependent trend in the data. For these reasons, we prefer to leave our results in 
terms of observational, rather than modeled, constraints. Using the data we provide, readers 
can straightforwardly test the effects of modeling different correction factors. 

Another topic that has received somewhat less attention is the possibility that the vari- 
ability signal could decrease with redshift because different spectral regions are being shifted 
across the Chandra bandpass. We have performed a basic study of how spectral variation 
may change as a function of photon energy. In our case ( §3.81 and §4.31) . we found some 
evidence that hard spectral regions may show increasing influence of a constant component 
that shifts at higher redshifts into an energy range where Chandra is most sensitive. At 
our current stage of understanding, attempting to correct for this effect would introduce 
systematic uncertainty, and we prefer to stay with our empirical approach. 

Increased variability at higher redshifts would have intriguing physical implications, al- 
though we note that it is not obv iously expected in the currently-favored "cosmic downsizing" 



pictu re for SMBH growth (e.g., ICowie et al.l 120031 : iMarconi et al.ll2004l : iBrandt fc Hasinger 



20051 . and references therein). In this picture, we observe multiple largely- independent gen- 
erations of growing SMBHs as we look back in cosmic time, with the more massive SMBHs 
generally growing earlier. At a fixed high luminosity, we then do not expect to see SMBH 
masses dropping (and thus perhaps variability strength rising) toward higher redshifts. 



4.3. Spectral Variation 

Quasar variation appears similar to that observed for lower-luminosity Seyfert AGN 
in several respects. The X-ray spectra of quasars vary in shape as well as brightness. In 
general, the spectrum of a given source is flatter when that source is in a fainter state, 
and the spectrum steepens as the source brightens. The change in HR is anti-correlated 
with the change in luminosity (Figure [T6|l . demonstrating the steepening-brightening trend. 
Direct spectral fits for multiple epochs of high-S/N sources (Figure [T7|) demonstrate a strong 
anti-correlation between AF and F that tends to move spectra toward F ~ 2. 

The power law spectral index (derived from our simple fit model) may vary significantly 
for a given source. Based on our fits to sources with large numbers of counts, intrinsic 
variation by AF ~ 0.1 — 0.2 is common. Some caution is therefore warranted when using 
the spectral slope as an indicator of quasar physical properties such as black hole mass or 
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accretion rate, as the spectral slope (measured according to our simple model) can vary 
with X-ray luminosity over time. Of course, the variation mechanism may contribute to a 
relationship between T and accretion rate {L/LEdd, where LEdd is the Eddington rate) that 



has been found in s ingle-epoch observations of quasar ensembles (e.g., IShemmer et al 



2008 



Risaliti et al.l 120091 ). For the epochs in Figure [T8| we measure a best-fit relationship: 



AF = (0.534 ± 0.104) log(iVi/iVo) - (0.039 ± 0.012), 



(10) 



where Ni/Nq is the ratio of power law normalizations between the two epochs, and AF = 
Fi — Fq is the change in F. (We note that this fit may not adequately describe the properties 
of individual sources, as we are simply assuming that one model describes the variation in all 
epochs. We do not have sufficient data to fit each source individuall y.) The slope in Equa- 
tion [TU] is a bit steeper than the slope (0.31 ±0.01) in Equation (1) of IShemmer et al.l (|2008[ ) 
that related F to L/ LEdd for individual AGN. The ove r all effect may be similar to the "intrin- 
sic" a nd "global" Baldwin effects (e.g.. lBaldwinlll977l : iKinney et al.lll990l : iPogge &: Peterson 
19921 ). in that the variation in an individual object may be more extreme than the overall 
trend observed across single-epoch measurements of a large number of sources. At the least, 
intrinsic variation contributes to scatter observed in the sample of sources. 

Simple absorption models did not adequately describe the observed anti-correlation 
between spectral hardness and luminosity, although we cannot rule out the possibility of more 
complex, multi-parameter absorption models. However, a two component emission model 
did reproduce the trend reasonably well without fine-tuning. The emission model consists 
of a variable power law with F = 2 and a hard (F = 1.3), constant component normalized 
to ~ 20% of the variable power law flux at 1 keV. Previously, Pao04 suggested that the 
increasing contribution of a flatter, reflected component may result in lower variability levels 
for sources with "absorbed," or harder, spectra. In our data, the relationship between hard- 
and soft-band X-ray fluxes, extrapolated to zero soft-band flux, also suggests that a constant, 
hard spectral component is present in our ensemble of sources. Similar mode ls have been 
used to describe flux variation in Seyfert AGN such as MCG -6-30-15 (e.g 
2OO3I : IVaughan fc FabiaJ2004h . 



Taylor et al. 



If a typical quasar spectrum includes a non-negligible constant component that is flatter 
than the dominant component with power law index F = 2, we might expect that this 
component changes the observed spectral shape for higher-redshift quasars. There has been 
disagreement over whether spectral sha pes evolve with r edshi ft, but so far strong evidence 
for evolution has not been discovered. IShemmer et al.l (120051 ) flnd that the spectral slope 
for a joint fit to 10 of the most luminous quasars at4<2;<6.3isF^2, similar to that 
of lower-redshift quasars. However, this result applies to the resi-frame 2-10 keV range, 
and not necessarily to a harder component that becomes more evident at higher energies. 
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Vignali et al.l ([20051) report a slightly harder (F = 1. 9^0:09) obse rved-irame photon index for 



a joint fit of 48 quasars detected in X-rays at 4.0 < z < 6.3. iGreen et al.l ( 120091 ) find no 
evidence for spectral evol u tion i n a large sample of optically-selected quasars out to z < 5.4. 
However, iBechtold et al.l ( l2003l ) previously obser ved spectral f latten ing in a sample of 17 
optically-selected quasars at 3.7 < z < 6.3, and iKelly et al.l (120091 ) observed marginally- 
significant spectral flattening with redshift out to ^ ~ 4.7. 

To determine how a two-component model would affect the observed spectrum at high 
redshift, we randomly generated spectra for 9 Chandra sources aX z = 4.2. We simulated 
each observation from a model with two power laws having F = 2 and F = 1.3. The flatter 
power law was normalized to 20% of the steeper component at rest-frame 1 keV; it starts to 
rise above the steeper component at rest-frame 10/(1 -|- z) keV. Then we fit each spectrum 
with a power law model to determine typical measured photon indices, finding an average 
value of F ^ 1.72 ± 0.17. 

For comparison, we fit Galactic-absorbed power law mode ls to higher-energy reg ions of 
spectra for quasars at 2; > 4. The quasars were described in IShemmer et al.l (120051 ). with 
spectra kindly provided to us by O. Shemmer. We fit the ohserved-haiae energy range £"0- 
Ei keV, where Eq = 10/(1 + z) keV and Ei varied, in order to characterize the region where 
a hard component may dominate (in our simple model). We fit the spectra of all three EPIC 
cameras simultaneously, requiring them to have the same photon index, but allowing slightly 
different normalizations to account for cross-calibration differences between instruments. In 
all cases, the photon indices for the observed frame spectra were flatter than for the rest- 
frame 2-10 keV region. For a fit range extending up to Ei = 5 keV in the observed frame, 
the spectrum flattened by AF = 0.1 to 0.9. For higher values of Ei, the spectrum appeared 
to flatten more, although background contamination becomes higher at these energies as 
well. This result suggests that a hard spectral component is becoming more evident at high 
energies. Although inspection indicates that backgrounds are generally well below count 
rates in our fit region, it would be useful to obtain larger, high-quality samples to examine 
the shape of the high-energy quasar spectrum more reliably. 

In addition to the hypothesis of a constant harder component in quasar spectra, the 
assumption (based on studies of local Seyfert AGN) that the relationship between soft and 
hard X-ray count rates can be extrapolated linearly should also be tested with larger, high- 
quality data sets. The two-power-law model we use here is the simplest attempt to explain 
the putative constant component, but more sophisticated models may be required that cut 
off at higher energies if the flatter component is not evident in high-redshift spectra. Our 
results also suggest that studies of high-redshift populations should test multi-component 
models against individual sources in their data to constrain spectral complexity beyond a 
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single power law shape. The spectral shape measured from a stack or joint fit can be biased 
toward the brighter sources in a sample, so care should ideally be taken to account for 
selection effects and brightening-steepening behavior as sources vary. 



4.4. [O III] Emission Strength 



In our sample, sources classified as X-ray non- variable have [O III] A5007A emission 
lines that are about 12A stronger (on average) than those of variable sources ( §3.10p . This 
was the most prominent effect we found when comparing composite spectra of variable 
and non-variable AGN. In the optical bandpass, variability has been observed to be related 
to both [O III] and H/j emission. Quasars with stronger [O III] emission have /arg'er optical 
variability amplitudes (IGiveon et al.lll999uMao et al.ll2009l ). although with signficant scatter. 
The cause of this relation is not known; it has been suggested that the effect may be due 
to an enhanced ionization rate from a variable ionizing continu um (IGiveon et al.lll999l ). or 
perhaps associated with accretion instabilities or star formation ( iMao et al.ll2009l ). However, 
in the X-rays, we see the opposite effect: stronger [O III] emission is associated with lower 
levels of X-ray variability. 

The effect may be anomalous. The set of sources for which we can sensitively measure 
[O III] emission and X-ray variability is relatively small, and the effect is not highly sig- 
nificant. However, we do not see a significant difference (according to a KS test) between 
the distributions of monochromatic luminosity {uL^ at 5100A), X-ray count-rate luminosity, 
or mean numbers of X-ray counts for the variable and non-variable sources. There may be 
other factors differentiating the two subsamples, or the effect may be a random deviation in 
our small sample. 

In our search for short term variability on kilosecond time scales ( §3.9p . we identified an 
unusual quasar, SDSS J1238, that was flagged as variable in multiple epochs. This quasar is 
NLSl-like, with atypically weak [O III] A5007A emission, a B.fsX48Q2 FWHM f«1860 km s-\ 
and stronger ionized Fe emission compared to the average SDSS quasar. A tende ncy for rapid 



1996; 


Pogge 


2000; 


Komossa 


2007) 



With EWs from a few to ^lOOA, the [O III ] A5007A lines in our study span the range of 
EWs originally used by lBoroson fc GreenI (119921 ) to define "Eigenvector 1." This empirically- 
based collection of properties includes an anticorrelation between Fe II and [O III] emission 
strengths; it also includes trends with radio-loudness and some H^ line features such as 
FWHM and asymmetry. NLSl AGN fall at one end of this eigenvector, generally having 
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narrow lines, weak [O III] emission, and strong ionized Fe emission. The quasar J1238 also 
shows similar properties when compared to typical SDSS quasars. In our sample. X-ray "non- 
variable" quasars may exemplify the (radio-quiet) population at the opposite (strong- [O III]) 
end of Eigenvector 1, while quasars showing moderate X-ray variability have [O III] emission 
levels more typical of ordinary SDSS quasars^ The sample of strong- [O III] emitters that 
have been observed multiple times in X-rays should be expanded in order to test the potential 
relation between X-ray variation and Eigenvector 1 more sensitively. Some new sources can 
be added from the archives as X-ray and optical surveys progress, but targeted observations 
would be most effective to obtain a significant ensemble of the strongest- [O III] emitters. 

A link between X-ray temporal properties and [O III] emission would be interesting 
because it would represent another physic al connection between small-scale (disk corona) 



and large-scale (NLR) AGN physics (e.g., iBrandt fc BoUerl Il998l ). (Of course, with larger 



data sets we could also test more sensitively whether additional parameters, such as Balmer 
line profiles, are related to X-ray variation.) In one category of models, a third parameter 
modulates both X-ray and NLR emission, creating an indirect relation between the two. For 
example, [O III] em ission strength is believed to be dominated by geometric factors (e.g.. 



Baskin fc Laorll2005l ). The Eddington ratio (L/LEdd) of bolometric to Eddington luminosity 



could control both the obscuration of X-ray/UV radiation essential for NLR ionization (e.g.. 



Abramowicz et al.lll980l : iBoroson fc Greenlll992l : IChang et al.ll2007l ) and also X-ray emission 



properties. Alternatively, some physical effect could create a more direct connection between 
the corona and NLR. Very large NLRs that produce the strongest [O III] emission lines 



interact with the NLRs of radio- quiet Seyfeit AGN (e.g. 



Leipski et al. 


2006; 


Wane et al. 


2011) 



Falcke et al. 


1998; 


Ho & Peng 


2001 



vector 1 properties, including radio loudness and X-ray spectral slope, and radio-loudness 
(indicative of jet activity) is associated with decreased variability in our subsample of radio- 
loud quasais ( §3.3.ip . 



^°0n the long timescales we gene rally sample, it i s not immediately clear how X-ray variability should 
relate to Eigenvector 1. For example. iMcHardv et all (|2007() find a break in the power spectral density (PSD) 
fmiction at w 10^^ Hz for the NLSl Ark 564, below which its PSD drops significantly. However, the black 
hole masses for our sample are generally much higher than those for which PSD studies have been performed 
(Figure [5]). We therefore expect any PSD breaks for quasar- luminosity NLSl analogs in our sample to be at 
lower frequencies of perhaps ~ 10~^ Hz, which corresponds to longer timescales than we adequately sample 
in this work. 



5. SUMMARY AND CONCLUSIONS 



In this study, we have examined archived Chandra X-ray observations of 264 SDSS 
spectroscopic quasars to search for X-ray variability, characterize it, and test whether this 
variability is related to other quasar properties. Our findings include the following: 

1. We find strong evidence of X-ray variation in ^30% of the quasars in our sample 
overall. Our sensitivity to variation increases with the number of source counts; 70% 
of sources with >1000 counts per epoch are detected as variable. 

2. Quasars in our sample typically vary with a standard deviation of fractional variation 
of ~16%. This amount of variation is not large enough to explain the scatter in 
X-ray-to-optical ratios as being due to variation alone. On time scales shorter than a 
few xlO^ s, the ensemble variability falls below our detection limit. Coupled with the 
flatter trend of variability on longer time scales, this suggests a "break" in the trend 
at ~ (2 - 5) X 10^ s. 

3. We find no evidence that higher- redshift quasars are more variable than lower-redshift 
quasars, as has been suggested in previous studies. However, this analysis is compli- 
cated by the fact that we do not have many quasars at high redshift with large numbers 
of counts (^100) to enable sensitive tests. 

4. The X-ray spectra of quasars tend to be flatter when fainter and steeper when brighter, 
as is seen in the case of some local Seyfert AGN. We were able to reproduce this trend 
with a simple, two-parameter power-law model that has been used to describe Seyfert 
variability. Spectral fits to bright sources show an anti-correlation between AF and F; 
quasar spectra tend to flatten or steepen as necessary to bring them back to F ~ 2 as 
they vary. 

5. As soft-band count rates are extrapolated to zero, a significant hard-band flux remains. 
This suggests that quasar spectra have an underlying constant, hard spectral compo- 
nent, following the model proposed for some Seyfert AGN. The constant fraction of 
the hard-band count rate (measured in observed-frame bands) likely increases with 
redshift as different segments of the constant and variable spectral components shift 
through the bandpass. 

6. A search for intra-observation variation on time scales of 1 ks revealed one unusual 
source, J1238, with strong, short-term variability. The optical spectrum of J1238 
shows it is an NLSl-type object. For the full sample, we constrain the rates of signif- 
icant variation in 1 ks bins to be < 37.3 yr~^ (absorption/dimming) and < 51.5 yr~^ 
(emission), using observed- frame years. 
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7. We generated upper limits on the rate of observations showing at least a magnitude 
F of fractional variation, where an "observation" is representative of the epochs in 
Sample HQ. As examples, \F\ > 100% is rare, occuring <4% of the time. \F\ > 25% 
occurs in fewer than 25% of observations. 

8. Median spectra suggest that sources with higher (detectable) levels of X-ray variability 
have weaker [O III] emission. Additional data are required to confirm the relation 
between Eigenvector 1 properties such as [O III] strength and X-ray variability, and 
(more generally) to test how this phenomenon could connect small-scale (corona) and 
large-scale (NLR) AGN structures. 



The sample of serendipitously-observed quasars continues to expand for variability stud- 
ies. The Chandra archive and spectroscopic quasar catalogs continue to grow over time. 
Incorporating ROSAT PSPC observations would extend the observed- frame time baseline 
to >20 yr, while archived XMM-Newton data can greatly increase sample sizes and provide 
simultaneous optical/UV monitoring. eROSITA, scheduled to launch within the next two 
years, will provide a sensitive new survey s canning the X-ray sky multiple times, with each 
scan taking about 6 months to complete (jCappelluti et al.l 120111). The SDS S-III Baryon 
Oscillation Spectroscopic Survey (BOSS) quasar survey (e.g.. lRoss et al.ll201ll ) is more than 
doubhng the current SDSS spectroscopic quasar sample, with most of the new quasars hav- 
ing redshifts z > 2.2. Although these sources will generally be fainter, they will provide 
much-needed leverage for studies of accretion evolution over cosmic time. In addition to the 
sorts of analyses we have conducted in this paper, these expanded data sets may permit 
us to place detailed constraints on the quasar variability power spectrum by comparing the 
data to light curves that are simulated from different PSD models. 

New AGN are increasingly being identified by temporal properties, permitting current 
and future surveys to go beyond the quasar realm to consider AGN that are blended with 
host galaxy emission. However, there is a great need for new data-mining and statistical 
techniques that will appropriately characterize the properties of fainter sources in these new 
surveys, while accounting for instrumental cross-calibration and perhaps selection biases in 
the samples. 
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Table 1. Properties Common to Each Source 



SDSS Name z Number Known^ Known ^2500a'' Const Rate Is Const Rate HQ Is 

(J2000) of Obs RL BAL IQ-^^ erg lO^^ cts Var 10"^ cts Var 

cm^^ s^^ Hz^^ cm^^ s^^ cm^^ s^^ HQ 
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- 000424.4 


1.038 


2 


1 





0.444 


219.874 + 16.620 
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2 







1.584 


3.963 ± 2.731 
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+ 002013.1 
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2 


* 





0.587 


3.226 + 1.189 





0.323 + 0.119 
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0.267 
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1 
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2 







0.700 
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024230.65 


- 000029.6 


2.505 


2 




1 


1.195 


0.597 + 0.352 


1 


0.060 + 0.035 


1 


074408.41 


+ 375841.1 


0.881 


4 







1.933 


3.429 + 1.572 


1 


0.343 + 0.157 


1 


074417.47 


+ 375317.2 


1.067 


5 


1 





2.541 


148.129 + 8.189 


1 


14.988 + 0.738 





074502.90 


+ 374947.0 


0.593 


4 


* 





0.349 


30.388 + 5.093 









074524.97 


+ 375436.7 


0.406 


4 


* 





0.671 


19.774 + 4.516 


1 






074545.01 


+ 392700.8 


1.629 


2 







0.875 


22.872 + 3.357 









075502.11 


+ 220346.8 


0.400 


2 







0.799 


9.354 + 1.006 





0.935 + 0.101 





080731.76 


+ 211754.3 


1.194 


2 







1.175 


19.689 + 4.555 


1 


1.969 + 0.455 


1 


080749.15 


+ 212122.3 


2.232 


2 







1.537 


10.819 + 3.168 





1.082 + 0.317 





081426.45 


+ 364713.5 


2.732 


2 


1 


1 


0.375 


4.104 + 2.895 





0.410 + 0.290 





083454.89 


+ 553421.1 


0.241 


2 


1 





0.733 


37.978 + 3.691 





3.798 + 0.369 





083633.54 


+ 553245.0 


1.614 


2 




1 


1.421 


18.285 + 6.333 









084207.58 


+ 322646.7 


1.685 


2 







6.344 


49.105 + 4.445 









084308.18 


+ 362439.6 


1.502 


2 


* 





0.412 


9.442 + 3.151 





0.944 + 0.315 





084905.07 


+ 445714.7 


1.259 


2 


* 





0.333 


19.417+1.425 





1.942 + 0.143 





084943.70 


+ 450024.2 


1.592 


2 







1.820 


30.711 + 2.160 









090900.43 


+ 105934.8 


0.162 


2 







2.984 


722.397 + 30.768 


1 


72.240 + 3.077 


1 
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090928.50 H 


h 541925.9 


3.760 


2 


* 


1 


0.768 


0.641 ± 0.732 









091029.03 H 


h 542719.0 


0.526 


2 







1.274 


156.138 ±4.701 


1 


15.614 ±0.470 


1 


091127.61 H 


h 055054.0 


2.793 


2 




1 


3.561 


21.430 ±2.804 









091210.34 H 


h 054742.0 


3.241 


2 







3.348 


5.948 ±3.310 









091752.54 H 


h 414530.5 


1.277 


2 







0.777 


8.723 ±3.373 









092108.62 H 


h 453857.3 


0.174 


2 


* 





1.476 


538.732 ± 21.725 


1 






092314.48 H 


h 510020.7 


1.388 


2 







1.323 


17.125 ± 10.509 









094745.14 H 


V 072520.6 


0.086 


2 







5.939 


320.078 ± 8.782 





32.008 ± 0.878 





095240.16 H 


h 515249.9 


0.553 


2 







1.511 


92.083 ± 9.037 





9.208 ± 0.904 





095243.04 H 


h 515121.0 


0.862 


2 







2.522 


110.763 ± 10.977 









095542.12 H 


h 411655.2 


3.420 


3 


* 





0.886 


3.037 ± 7.793 


1 


0.302 ±0.215 


1 


095544.91 H 


h 410755.0 


1.921 


3 


* 


1 


0.710 


1.604 ±5.929 





0.160 ±0.104 





095548.13 H 


V 410955.3 


2.308 


2 







1.232 


17.079 ± 7.844 









095640.38 H 


h 411043.5 


1.887 


3 


* 





0.244 


2.368 ±6.652 


1 


0.237 ±0.181 


1 


095820.44 H 


V 020303.9 


1.356 


2 


* 





0.409 


2.876 ± 0.963 





0.288 ± 0.096 





095835.98 H 


h 015157.0 


2.934 


2 


* 







1.892 ± 0.965 





0.189 ±0.097 





095857.34 H 


h 021314.5 


1.024 


2 


* 





0.299 


91.785 ± 4.438 


1 






095858.68 H 


h 020138.9 


2.454 


4 







1.154 


19.584 ±3.237 


1 


1.958 ±0.324 


1 


095902.76 H 


h 021906.3 


0.345 


2 







0.671 


108.981 ± 4.718 





10.898 ± 0.472 





095918.70 H 


h 020951.4 


1.157 


4 


1 





0.426 


83.168 ±6.029 


1 


8.317 ±0.603 


1 


095924.46 H 


V 015954.3 


1.236 


4 







1.100 


53.427 ±5.123 





5.343 ±0.512 





095949.40 H 


h 020140.9 


1.753 


3 


* 





0.730 


19.564 ± 2.730 





1.956 ±0.273 





095957.97 H 


h 014327.3 


1.618 


2 


* 





0.399 


7.425 ± 1.990 









100012.91 H 


V 023522.8 


0.699 


5 







0.788 


31.959 ±6.387 


1 


2.948 ± 0.595 


1 


100014.13 H 


h 020054.4 


2.497 


5 


* 





0.588 


7.728 ±2.699 


1 


0.795 ±0.223 


1 


100024.39 H 


h 015053.9 


1.664 


6 


* 





0.510 


3.951 ± 2.383 


1 


0.400 ±0.161 





100024.64 H 


V 023149.0 


1.318 


5 







0.886 


21.383 ±5.368 


1 


2.138 ±0.537 


1 


100025.24 H 


h 015852.0 


0.373 


4 







1.060 


141.057 ±8.725 


1 


14.106 ±0.873 


1 


100043.13 H 


h 020637.2 


0.360 


6 


* 





0.632 


19.413 ±5.233 


1 


1.941 ±0.523 


1 


100055.39 H 


h 023441.3 


1.403 


4 


* 





0.514 


15.301 ± 3.707 





1.559 ±0.322 





100058.84 H 


h 015400.2 


1.559 


5 


* 





0.386 


18.397 ±5.385 





1.840 ±0.538 





100104.32 H 


h 553521.6 


1.537 


2 







1.247 


16.386 ± 3.820 





1.639 ±0.382 





100114.29 H 


h 022356.8 


1.799 


4 







0.876 


13.170 ±2.702 


1 


1.317 ±0.270 


1 


100116.78 H 


h 014053.5 


2.055 


3 


* 





0.529 


10.119 ±2.994 





1.182 ±0.268 





100120.26 H 


h 023341.3 


1.834 


2 


* 





0.365 


4.476 ± 1.415 





0.448 ±0.142 





100130.37 H 


h 014304.3 


1.571 


5 


* 





0.312 


4.226 ±3.122 





0.407 ±0.277 


1 


100145.15 H 


h 022456.9 


2.032 


2 


* 





0.311 


4.832 ± 1.187 





0.483 ±0.119 





100201.51 H 


h 020329.4 


2.008 


6 







1.232 


1.998 ±4.013 





0.210 ±0.247 





100205.36 H 


h 554257.9 


1.151 


2 







1.834 


16.737 ± 5.599 









100219.49 H 


h 015537.0 


1.509 


2 


* 





0.428 


17.506 ± 3.740 









102350.94 H 


V 041542.0 


1.809 


2 


* 





0.766 


1.854 ± 1.639 


1 


0.185 ±0.164 


1 


103222.85 H 


h 575551.1 


1.243 


2 


1 





0.330 


22.592 ± 2.706 









103227.93 H 


h 573822.5 


1.968 


3 


* 





0.287 


34.075 ± 3.835 





3.407 ±0.383 





104829.95 H 


h 123428.0 


0.442 


3 







0.790 


49.208 ±4.251 


1 


4.921 ± 0.425 


1 
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105015.58 H 


h 570255.7 


3.273 


2 


* 





0.450 


8.182 ± 7.308 









105039.54 H 


V 572336.6 


1.447 


2 


1 





0.694 


23.350 ± 8.775 









105050.14 H 


h 573820.0 


1.281 


2 







0.797 


26.146 ±8.842 









105239.60 H 


h 572431.4 


1.112 


3 







2.532 


51.458 ± 12.947 


1 






105316.75 H 


h 573550.8 


1.205 


3 







0.684 


81.026 ± 15.179 





8.103 ± 1.518 





105518.08 H 


h 570423.5 


0.696 


2 







1.123 


47.098 ± 12.371 





4.710 ± 1.237 





111354.66 H 


h 124439.0 


0.680 


2 


1 





0.107 


30.643 ± 19.016 


1 






111422.47 H 


h 531913.2 


0.885 


3 







1.187 


38.261 ± 10.991 


1 


3.826 ± 1.099 


1 


111452.84 H 


h 531531.7 


1.213 


3 







0.688 


43.474 ± 11.329 





4.347 ± 1.133 





111518.58 H 


h 531452.7 


1.540 


3 







1.735 


16.952 ±5.661 





1.695 ± 0.566 





111520.73 H 


h 530922.1 


0.877 


3 







2.082 


0.183 ±2.450 





0.018 ±0.245 





111816.95 H 


h 074558.1 


1.736 


2 







13.013 


114.067 ±6.977 


1 






111840.56 H 


h 075324.1 


1.463 


2 







0.888 


24.815 ±4.632 









111946.94 H 


h 133759.2 


2.023 


2 







1.127 


8.410 ± 2.999 









112026.20 H 


V 134024.6 


0.982 


2 







0.950 


68.886 ± 11.089 


1 






112048.99 H 


h 133821.9 


0.513 


2 







0.793 


3.115 ± 1.099 





0.312 ± 0.110 





112106.07 H 


h 133824.9 


1.944 


2 







1.152 


45.943 ± 14.019 









112213.65 H 


V 041548.7 


3.517 


4 


* 





0.406 


1.305 ±4.408 





0.142 ± 0.321 





112320.73 H 


V 013747.4 


1.469 


2 







17.207 


108.419 ± 10.153 









112404.52 H 


h 040418.1 


3.841 


5 


* 





0.929 


2.392 ± 5.626 





0.228 ±0.459 





114636.88 H 


h 472313.3 


1.895 


2 


1 





1.366 


51.248 ±4.892 





5.125 ± 0.489 





114651.21 H 


h 471732.5 


3.130 


2 


* 


1 


0.692 


0.999 ± 2.414 









114656.73 H 


h 472755.6 


0.668 


3 







0.872 


109.369 ±8.629 





10.937 ± 0.863 





115324.46 H 


V 493108.7 


0.334 


2 


1 





2.854 


1481.928 ± 22.813 


1 






115838.56 H 


h 435505.8 


1.208 


2 


1 





0.434 


15.222 ±3.356 





1.522 ± 0.336 





115911.43 H 


V 440818.3 


1.438 


2 







1.320 


22.118 ±6.177 









120104.66 H 


V 575846.9 


1.842 


2 







2.019 


26.017 ±8.028 


1 


2.602 ± 0.803 


1 


120106.14 H 


h 580336.6 


1.087 


2 


* 





0.470 


20.459 ±6.911 





2.046 ± 0.691 





120233.39 H 


V 580501.8 


3.424 


2 


* 





1.389 


11.738 ± 4.683 





1.174 ±0.468 





120924.07 H 


h 103612.0 


0.395 


5 




1 


6.871 


7.827 ± 26.396 









120924.80 H 


V 102553.9 


0.263 


2 


* 





1.072 


110.776 ± 38.543 









120937.02 H 


V 103756.9 


1.994 


5 







1.102 


6.221 ± 23.135 









120949.46 H 


h 102146.8 


2.312 


6 


* 





0.898 


12.216 ±35.370 









120959.05 H 


V 104320.1 


1.314 


4 







0.855 


13.081 ± 404.062 









121013.57 H 


h 104853.7 


1.080 


4 


1 





0.912 


52.349 ±35.155 









121111.46 H 


h 100826.8 


1.994 


2 







0.952 


15.664 ±9.156 









121342.95 H 


h 025248.9 


0.641 


2 







0.674 


14.385 ±3.101 





1.438 ±0.310 





121440.27 H 


h 142859.1 


1.625 


2 




1 


3.824 


1.143 ±3.444 









122418.00 H 


h 070949.2 


0.981 


7 







0.829 


5.848 ±4.537 


1 


0.497 ±0.211 


1 


122448.15 H 


h 125413.3 


1.062 


2 







2.384 


0.639 ± 0.454 





0.064 ± 0.045 





122511.91 H 


h 125153.6 


1.255 


3 







2.188 


86.787 ± 4.564 


1 


8.679 ± 0.456 


1 


122515.65 H 


h 124441.0 


1.664 


2 







1.046 


18.494 ± 2.974 


1 


1.849 ± 0.297 


1 


122722.12 H 


h 075555.0 


3.168 


2 


1 





0.986 


8.822 ±8.352 









122826.33 H 


h 130106.2 


3.229 


3 


* 





0.816 


3.757 ± 5.737 
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122923.73 


+ 075359.2 


0.854 


2 







0.838 


70.518 ± 6.482 


1 


7.052 ± 0.648 


1 


123320.92 


+ 110702.4 


1.206 


2 







1.107 


15.059 ± 7.321 









123346.21 


+ 130905.7 


1.368 


2 







0.746 


17.454 ± 6.481 









123410.72 


+ 111732.6 


0.817 


2 







1.556 


39.691 + 8.514 


1 


3.969 + 0.851 


1 


123527.75 


+ 121338.8 


0.726 


4 







3.252 


71.251 ± 17.906 


1 


7.565 + 1.375 


1 


123540.19 


+ 123620.7 


3.208 


2 


* 





0.552 


2.839 ± 5.799 









123618.94 


+ 121010.0 


0.993 


2 







0.803 


5.519 + 6.425 









123622.94 


+ 621526.6 


2.587 


20 


* 





0.277 


7.122 + 3.914 


1 


0.712 + 0.391 


1 


123715.99 


+ 620323.3 


2.068 


6 


* 





0.511 


1.963 + 2.437 


1 


0.196 + 0.244 


1 


123759.56 


+ 621102.3 


0.909 


10 







1.309 


59.321 +8.443 


1 


5.932 + 0.844 


1 


123800.91 


+ 621336.0 


0.440 


15 







0.858 


36.941 +8.695 


1 


3.694 + 0.870 


1 


124107.11 


+ 113701.7 


1.413 


2 







1.093 


26.428 + 7.978 









124210.41 


+ 115223.8 


0.298 


2 







0.689 


72.599 + 8.249 


1 






124255.31 


+ 024956.9 


1.459 


4 







0.770 


22.482 + 4.444 


1 


2.248 + 0.444 


1 


125849.83 


- 014303.3 


0.967 


8 







5.162 


136.037 + 16.400 


1 


13.604+ 1.640 


1 


125919.26 


+ 124829.0 


0.701 


2 







0.839 


26.000 + 4.759 


1 






130216.13 


+ 003032.1 


4.468 


2 


* 





0.618 


0.498 + 0.912 









132852.11 


+ 472218.3 


1.932 


2 







0.874 


8.923 + 3.570 









132938.57 


+ 471854.6 


1.027 


2 


* 





0.316 


26.418 + 10.304 









133004.72 


+ 472301.0 


2.825 


3 


1 


1 


0.805 


1.095 + 5.995 









133223.26 


+ 503431.3 


3.807 


2 







1.922 


6.185 + 2.478 





0.618 + 0.248 





134425.94 


- 000056.2 


1.096 


2 







1.362 


5.921 + 4.200 









135854.44 


+ 623913.1 


1.228 


3 







1.523 


63.993 + 6.630 


1 






140041.11 


+ 622516.2 


1.878 


2 


* 





0.586 


9.450 + 3.582 





0.945 + 0.358 





140146.53 


+ 024434.7 


4.441 


2 


* 





1.957 


9.610 + 2.363 


1 






140354.57 


+ 543246.8 


3.258 


2 


* 





0.508 


4.817 + 1.396 


1 






141500.38 


+ 520658.5 


0.424 


6 







1.016 


48.668 + 9.485 


1 


4.867 + 0.949 


1 


141533.89 


+ 520558.0 


0.986 


6 







1.178 


28.762 + 6.984 


1 


2.876 + 0.698 


1 


141551.27 


+ 522740.6 


2.585 


2 


* 





0.503 


3.131+2.255 





0.313 + 0.226 





141551.59 


+ 520025.6 


1.514 


3 







0.980 


8.662 + 4.682 


1 


0.866 + 0.468 


1 


141642.42 


+ 521812.7 


1.284 


9 









18.385 + 7.574 


1 


1.839 + 0.757 


1 


141647.20 


+ 521115.2 


2.153 


5 







1.530 


16.566 + 5.502 





1.635 + 0.327 





141905.17 


+ 522527.7 


1.606 


3 


* 





0.447 


13.304 + 2.754 









141908.18 


+ 062834.8 


1.437 


2 


1 





7.273 


1305.172 + 45.235 


1 






142005.59 


+ 530036.7 


1.647 


17 


* 





0.389 


11.174 + 7.145 


1 


1.117 + 0.714 


1 


142015.64 


+ 523718.8 


1.674 


6 


* 





0.781 


14.832 + 6.621 









142052.43 


+ 525622.4 


0.676 


24 







1.969 


83.371 + 22.501 


1 


8.358 + 2.188 


1 


142147.09 


+ 532405.7 


3.039 


5 


* 





0.421 


3.556 + 2.422 


1 


0.356 + 0.242 


1 


142301.08 


+ 533311.8 


1.863 


2 







2.837 


22.557 + 2.538 





2.256 + 0.254 





142305.04 


+ 240507.8 


4.105 


2 


* 


1 


0.561 


4.479 + 1.408 





0.448 + 0.141 





142455.69 


+ 351356.6 


1.255 


2 







0.861 


5.769 + 3.806 









142507.32 


+ 323137.4 


0.478 


2 


1 





1.544 


79.144 + 14.726 









142530.09 


+ 335217.3 


1.185 


2 







0.768 


5.538 + 4.512 





0.554 + 0.451 





142532.83 


+ 330124.9 


1.200 


2 


1 





0.286 


22.939 + 8.534 





2.294 + 0.853 
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142539.01 H 


h 331009.4 


2.306 


2 







1.881 


35.656 ± 9.731 





3.566 ± 0.973 





142543.30 H 


h 335543.6 


1.133 


3 







0.675 


38.937 ± 11.438 





3.894 ±1.144 





142545.53 H 


h 332603.3 


2.963 


2 


* 





0.469 


3.237 ± 4.394 





0.324 ±0.439 





142551.17 H 


h 350113.0 


0.898 


2 







1.164 


28.423 ± 8.816 





2.842 ±0.882 





142557.63 H 


h 334626.2 


0.351 


2 







2.016 


23.779 ± 8.883 





2.378 ±0.888 





142620.30 H 


h 351712.1 


1.748 


2 




1 


1.364 


2.221 ± 6.837 


1 


0.222 ±0.684 


1 


142622.66 H 


V 334202.3 


1.349 


2 







2.069 


35.988 ± 11.628 









142623.15 H 


h 351154.9 


3.503 


2 


* 





0.561 


0.879 ±3.865 





0.088 ±0.387 





142640.83 H 


h 332158.7 


1.542 


2 




1 


0.891 


2.121 ± 5.310 





0.212 ±0.531 





142730.19 H 


h 324106.4 


1.776 


2 







0.955 


1.899 ±3.277 





0.190 ±0.328 





142734.80 H 


h 352543.4 


0.340 


2 







1.303 


37.178 ± 10.303 





3.718 ± 1.030 





142738.21 H 


h 351132.1 


1.208 


3 







1.389 


15.525 ± 9.494 





1.664 ±0.726 





142738.36 H 


V 325320.0 


0.822 


2 


* 





0.574 


89.634 ± 10.516 





8.963 ± 1.052 





142810.31 H 


h 353847.0 


0.804 


2 







1.006 


99.180 ± 10.800 





9.918 ± 1.080 





142813.94 H 


V 334759.6 


2.243 


3 







0.967 


1.590 ±5.479 


1 


0.256 ±0.509 


1 


142817.81 H 


h 354021.9 


0.335 


2 







0.791 


59.363 ± 10.160 





5.936 ± 1.016 





142848.32 H 


h 350315.5 


2.111 


3 







1.049 


5.417 ± 5.954 





0.542 ±0.595 





142858.01 H 


y 344149.9 


3.076 


2 


* 





0.505 


2.851 ±4.813 





0.285 ±0.481 





142910.22 H 


h 352946.8 


2.224 


3 


1 





1.181 


9.653 ±4.223 





0.965 ±0.422 





142911.17 H 


h 330941.3 


1.116 


2 







0.779 


14.136 ±6.838 





1.414 ±0.684 





142912.87 H 


h 340959.0 


2.229 


4 







2.063 


14.479 ± 10.076 





1.468 ±0.895 





142915.19 H 


h 343820.3 


2.351 


2 







1.495 


10.811 ± 6.959 





1.081 ±0.696 





142917.20 H 


h 342130.3 


1.275 


2 







0.800 


60.240 ± 13.005 





6.024 ±1.300 





142942.64 H 


V 335654.7 


1.121 


3 


1 





0.961 


84.645 ± 15.627 





8.464 ±1.563 





142949.65 H 


h 324653.9 


2.175 


2 


* 





0.852 


8.675 ± 7.314 





0.867 ±0.731 





142954.70 H 


h 330134.7 


2.076 


2 







2.103 


20.133 ±8.187 





2.013 ±0.819 





143031.78 H 


V 330042.5 


1.071 


2 







1.212 


13.942 ± 8.611 





1.394 ±0.861 





143034.83 H 


h 335945.3 


1.115 


2 







2.628 


43.009 ± 9.442 





4.301 ± 0.944 





143106.77 H 


h 340910.8 


1.098 


3 







0.669 


20.320 ±9.116 





2.032 ±0.912 





143107.51 H 


h 342730.9 


4.270 


2 


* 





0.548 


2.269 ±4.713 









143132.13 H 


h 341417.3 


1.040 


2 







1.330 


57.982 ± 12.491 





5.798 ± 1.249 





143157.94 H 


h 341650.2 


0.715 


2 







4.634 


272.680 ± 22.804 





27.268 ±2.280 





143201.74 H 


h 343526.2 


1.071 


2 







0.934 


13.472 ± 8.661 









143219.53 H 


V 341728.9 


0.629 


2 







1.206 


106.076 ± 19.235 


1 


10.608 ± 1.924 


1 


143243.92 H 


h 330746.6 


2.088 


2 







1.621 


6.423 ±5.359 





0.642 ± 0.536 





143244.26 H 


h 350100.4 


1.038 


2 







0.957 


14.640 ± 7.434 


1 


1.464 ±0.743 


1 


143307.88 H 


V 342315.9 


1.950 


2 


* 





0.839 


19.433 ± 8.290 





1.943 ±0.829 





143331.79 H 


h 341532.7 


0.957 


2 







0.887 


1.822 ± 4.439 





0.182 ±0.444 





143335.68 H 


h 350133.1 


0.618 


2 







1.213 


3.466 ±5.645 









143345.10 H 


V 345939.9 


0.815 


2 


* 





0.472 


2.857 ±5.754 









143506.45 H 


h 335526.0 


3.940 


2 


* 





1.055 


6.056 ±4.582 





0.606 ± 0.458 





143547.62 H 


V 335309.6 


2.111 


3 


1 





0.315 


4.438 ± 6.344 





0.586 ±0.574 





143559.19 H 


h 334640.1 


0.948 


2 







0.899 


70.576 ± 14.211 





7.058 ± 1.421 





143604.64 H 


h 350428.5 


3.033 


2 


* 


1 


0.414 


0.033 ±2.619 





0.003 ± 0.262 
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Table 1 — Continued 



SDSS Name z Number Known^ Known ^2500a'' Const Rate Is Const Rate HQ Is 

(J2000) of Obs RL BAL 10-^7 erg 10"^ cts Var 10"^ cts Var 

cm~-^ s~^ Hz^^ cm~^ s~^ cm~^ s~^ HQ 



143617.81 


+ 353726.1 


1.448 


3 







1.016 


27.628 ± 12.945 





2.457 + 0.964 


1 


143624.30 


+ 353709.4 


0.767 


2 







1.151 


34.920 + 10.075 





3.492 + 1.007 





143624.61 


+ 352537.2 


1.060 


2 







1.052 


31.300 + 9.896 





3.130 + 0.990 





143626.63 


+ 350029.7 


1.242 


3 







0.915 


0.347 + 3.547 





0.035 + 0.355 





143627.78 


+ 343416.8 


1.883 


2 







1.026 


11.659 ± 6.829 





1.166 + 0.683 





143628.08 


+ 335524.3 


0.903 


3 







1.128 


30.806 ± 12.188 


1 


2.700 + 0.834 





143632.99 


+ 344253.4 


1.948 


2 







0.874 


14.846 ± 6.737 





1.485 + 0.674 





143651.51 


+ 343602.4 


0.296 


2 







0.877 


51.688 ± 13.078 


1 


5.169 + 1.308 


1 


143706.20 


+ 343659.2 


4.369 


2 


* 





0.600 


3.378 ± 6.450 









143841.95 


+ 034110.3 


1.737 


2 







2.346 


22.186 + 5.090 





2.219 + 0.509 





143859.05 


+ 033547.4 


0.734 


2 







1.291 


35.976 + 6.427 









144642.92 


+ 012552.4 


1.421 


2 







1.449 


14.527 + 6.972 









145206.45 


+ 580625.9 


1.440 


2 


* 





0.504 


36.140 + 2.998 


1 


3.614 + 0.300 


1 


145207.32 


+ 580454.7 


1.920 


2 


* 





0.301 


34.166 + 2.700 





3.417 + 0.270 





145215.59 


+ 430448.7 


0.296 


3 







0.927 


44.665 + 7.644 


1 


4.466 + 0.764 


1 


150407.51 


- 024816.5 


0.217 


2 


1 





2.657 


92.683 + 5.856 


1 






150948.65 


+ 333626.7 


0.512 


2 







0.568 


44.471 + 6.827 


1 


4.447 + 0.683 


1 


151413.58 


+ 553500.7 


1.319 


2 







1.584 


14.912 + 2.713 


1 






151451.28 


+ 552602.3 


1.842 


2 


* 





0.263 


6.818 + 1.632 





0.682 + 0.163 





151545.08 


+ 553518.4 


1.652 


2 


* 





0.412 


9.899 + 1.551 





0.990 + 0.155 





153308.65 


+ 301820.7 


4.455 


2 


* 





0.874 


-0.038 + 0.992 





-0.004 + 0.099 





155633.77 


+ 351757.3 


1.495 


2 


1 


1 


1.708 


21.838 + 3.413 









160410.22 


+ 432614.7 


1.538 


2 







2.301 


8.437 + 2.099 





0.844 + 0.210 





160630.60 


+ 542007.5 


0.820 


2 







1.085 


58.018 + 14.034 









160856.78 


+ 540313.7 


1.915 


2 


* 





0.767 


10.137 + 6.277 





1.014 + 0.628 





164025.02 


+ 464449.0 


0.537 


2 







1.604 


35.889 + 8.166 


1 


3.589 + 0.817 


1 


164733.23 


+ 350541.5 


0.861 


2 


* 





0.395 


38.988 + 8.173 





3.899 + 0.817 





165108.85 


+ 345633.7 


1.541 


2 







1.590 


8.492 + 3.861 





0.849 + 0.386 





170224.52 


+ 340539.0 


2.038 


2 







1.232 


5.852 + 3.762 





0.585 + 0.376 





170441.37 


+ 604430.5 


0.372 


2 







22.112 


402.461 + 12.583 


1 






171957.81 


+ 263416.6 


3.160 


2 


1 





0.399 


1.537 + 3.592 









172026.47 


+ 263816.0 


1.141 


4 







0.880 


17.377 + 44.859 





1.739 + 0.468 





172211.65 


+ 575652.0 


1.610 


2 


* 





0.342 


7.416 + 4.865 









173744.88 


+ 582829.6 


4.918 


2 


* 





1.041 


-0.130 + 3.845 









173801.22 


+ 583012.1 


0.330 


2 


* 





0.630 


73.887 + 12.107 





7.389 + 1.211 





173836.16 


+ 583748.5 


1.279 


2 







2.726 


4.692 + 4.609 









221453.84 


+ 140022.2 


1.523 


2 


* 





1.700 


21.435 + 3.355 





2.144 + 0.335 





221458.45 


+ 135344.7 


3.673 


4 


* 





0.672 


8.918 + 3.335 


1 


0.892 + 0.334 


1 


221738.41 


+ 001206.6 


1.121 


2 


* 





0.275 


6.937+ 1.132 


1 


0.694 + 0.113 


1 


221751.29 


+ 001146.4 


1.491 


5 


* 





0.450 


8.723 + 1.626 


1 


0.872 + 0.163 


1 


221755.20 


+ 001512.3 


2.092 


4 


* 





0.442 


8.789 + 1.548 


1 


0.879 + 0.155 


1 


232007.52 


+ 002944.3 


0.942 


4 







0.987 


87.608 + 9.669 





8.761 +0.967 





233130.08 


+ 001631.6 


2.659 


2 


* 





0.573 


4.565 + 2.527 





0.457 + 0.253 





235653.87 


- 010731.5 


0.601 


2 


* 





0.263 


35.650 + 5.004 





3.565 + 0.500 
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Table 1 — Continued 



SDSS Name 
(J2000) 


z Number 
of Obs 


Known'' 
RL 


Known 
BAL 


F 

2500A 
10-27 erg 


Const Rate 
10-" cts 

9 1 

cm s 


Is 
Var 


Const Rate HQ 
IQ-S cts 

— 9 — 1 

cm s 


Is 
Var 
HQ 





"Sources marked with an asterisk are not known to be radio-loud, but limits were not sensitive enough to guarantee 
that they were radio-quiet with high confidence. 

''Blank entries indicate where F2500 could not be reliably measured due to bad spectral bins. 



Table 2. Observations of Each Source 



SDSS Name 




Obsid 


Exposure 


OfE-Axis Angle 


In HQ 


TSTART^' 


Source Count Rate 


Count Luminosity 


Total Counts'' 


(J2000) 






(sec) 


(arcmin) 


Flag 


(sec) 


10"*^ cts cm-2 s-i 


10^2 cts s-i 


Full/Soft/Hard Band 


000622.60 - 000424, 


,4 


4096 


4451.707 


13.115 





176198559.771 


,g2+i4-7 
^^^^-14, 3 




274.582ll?;0f, / 192.9161^5:^29 / 85.2081^^3!^ 


000622.60 - 000424, 


,4 


5617 


16931.652 


7.950 


1 


238942682.443 


99Q+8.25 
^^^-8.02 


231.9301^-3^5 


960.227tlt f46 1 629.275l^»;237 / 343.528l?i:»^» 


004054.65 - 091526, 


,7 


4885 


9636.413 


7.774 


1 


210608981.312 


Q 7Q+2.17 
^•'■^-2.03 


294.766l233;^9| 


6.895l*:i35 / 5.81211^^° / 1.15llg:«^^, 


004054.65 - 091526, 


,7 


904 


38419.192 


14.503 





83056012.160 




519.462ll«^:f,^ 


119.043l}?-««i / 58.772l9-f3? / 63.200l98 ™ 


011513.16 + 002013, 


,1 


3203 


40581.117 


6.684 


1 


126613529.277 


9 in+0.782 
^■^='-0.772 




33.976+^*^3 1 24.072+«-*^/ 10.361+*-;^;^^ 


011513.16 + 002013, 


,1 


3204 


37616.425 


3.727 


1 


152506743.623 


4 62+°-^*^ 




42.287l?;gf3 / 35.5711^5 / 6.93ll*;^5,6 


014219.01 + 132746, 


,5 


1633 


1914.645 


14.098 





97032805.579 


fi 04+8-83 
D-U"l_8.45 


0.16llg:i« 


4.40ll3:«g/ 0.000lg:0^«/ 4.602l3:g8 


014219.01 + 132746, 


,5 


4010 


5064.646 


14.653 





161474857.473 


^■■^"-5.37 




10.006l*:|^3 / 3.427l?:3^5 / 6.912l*:^« 


014320.96 + 132429, 


,7 


1633 


1914.645 


12.922 





97032805.579 


n Q+6.83 
'-'■='-7.76 


-I7.96OI33I-,- 


4.69313:^?° / 2.579l?:*,^« / 2.2931?;°^^ 


014320.96 + 132429, 


,7 


4010 


5064.646 


14.251 





161474857.473 




A Sf.;zl+22.287 
'±.00'i_24 Qgg 


19.981™/ 8.035l^:3«/12.64ll53:?«« 



''The Chandra TSTART parameter indicates the time of the start of the observation in seconds since 1 Jan 1998. 

''The number of full-band counts is not exactly equal to the combination of soft- and hard-band counts due to band-dependent factors in aperture corrections and background estimation. 



Table 3. Observations of Each Source 



SDSS Name 


Obsid 


Source Counts^ 


BG Counts'' 


HR 


Approx. Lo.5-8'' 


(J2000) 




Full/Soft/Hard Band 


Full/Soft/Hard Band 




10**^ erg 



000622, 


,60 


- 000424, 


,4 


4096 


000622, 


,60 


- 000424, 


,4 


5617 


004054, 


,65 


- 091526, 


,7 


4885 


004054, 


,65 


- 091526, 


,7 


904 


011513, 


,16 


+ 002013, 


1 


3203 


011513, 


,16 


+ 002013, 


1 


3204 


014219, 


,01 


+ 132746, 


,5 


1633 


014219, 


,01 


+ 132746, 


,5 


4010 


014320, 


,96 


+ 132429, 


,7 


1633 


014320, 


,96 


+ 132429, 


,7 


4010 



254.269^ 



3+19,407 / 1 0R 989+16.471 / 7, p.r,o+11.461 
^-18.883 / -'-*°-^°^-15.695 / '-'-■'^^■^-11.251 

952.258tr3:f3? / 625.993lf,:f,g / 338.754l^Ji^^ 

oc:9+4.245 / c- „o9+3.996 / f, 070+2.730 
o.ooz_3 ggg / j.uoz_3 g24 / "■'-"°-2,754 

41.77111^:^^3 / 24.029li?:f,l / 18.773111;^°^ 
22.319+_lilt I 16.675i«;?l* / 5.94lll°?« 



40.886t?-^^^ 

9 CI 9+3,674 

1.385+5-219 
-2.480 



7.786 
6.972 
+2.171 
2.460 

^19 / 1 7(;c+3.509 / 
5.473 / ^- '"°-3.427 / 
+3.548 / 
3.579 



/ 34.795: 
/ -0.941^ 



/ 6.307; 



+4.338 / , 
4.930 / 1-1*° 



1.555t^;ig9 / 4.407lt:5M / -5-980 



+4.182 
3.725 

/ 3.61513-723 
-0.342t4;g|f 

■3.^'^_4.297 
+6.153 
7.003 



20.31311- / 6 
7.969l^iii / 

■2.856 



1.543 
77.272 



0,964 / 
10.485 / 
-7,605 / 
11 fic7+4.900 
li.D0(_2 ggy 

1.401+2-830 

1 



-0.933 
Q+2.799 
^-1.126 

8.62ll^-i^^ 

7 170+4.375 
(.1(0_2 214 



+4.287 / ,0 
-2.152 / 
+3.374 , , 
1.542 / * 
+2.591 / 
-0.608 / ^ 
-7.515 / 
-5.121 / 
+4.226 
2.377 
+2.606 
0.623 



633 
3.282 
0.780 

743 
7.397; 
0.776 

0.94llg;«? / 

\.m\t\f^^ I 7, 

1.392 



+3.147 
, -0.914 , 

21.536l|;gg2 / 3-628li:6i5 / 18.621 



.685 
.774 
.773 
.427 
/ 4.420 
/ 0.623 
.987 
254 
/ 5.767 





,363 






,833 




+3, 


,714 






,716 




.+2, 


,567 




'-0 


,548 




.+8, 


,395 






,362 




,+i. 


,620 
,654 


n c:q9+0.287 
u.ooz_Q 2ig 


,+2, 
'-0 


,491 
,461 


n 70/1+0.106 
-0-'84„o 


.+2, 


,580 






.695 




+4, 


,228 






,124 




.+4, 


,020 






.787 




+5, 


,903 






,459 





3216.715l23l 
3819.410li34- 
4854.19513850 
8554.465l|°24. 
298.165lio|; 
630.20lliig 

i+l 

s+517: 

-587, 



2.6543 
0.598j 
-295. 758"* 



-79.94311^ 



514 
886 
927 
112 
118 
038 
906 
192 
665 

334 rn 

666 y; 

718 
881 
716 
252 
361 
292 
938 
018 
323 



''The number of full-band counts is not exactly equal to the combination of soft- and hard-band counts due to band-dependent factors in aperture corrections and 
background estimation. 

''Calculated using a power law model with photon index F = 2. 
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Table 4. Simple Spectral Model for Bright Sources 



SDSS Name 


Obsid 


Norm 




r 


Edge Depth (r) 


(J2000) 




10 pn key cm s at 1 keV 










094745.14 + 072520.6 


7265 


6.26+°-«5 

— 0.318 


0, 


,091+0-0« 

—0.047 


3, 


4qi +0.469 
■^'-'^-0.416 


094745.14 


+ 072520.6 


6842 


6.37+°-349 

— 0.258 


0, 


,087+0-037 

—0.037 


3, 


146+0.357 
■ "-0.322 


095240.16 


+ 515249.9 


3195 


4.95+°-?f6 

"^" — 0.118 


2, 


00^+0.038 
"-*'-*'-^— 0.038 


0, 


■ooolg:«3 


095240.16 


+ 515249.9 


7706 


4.63+3,07 

— 0.654 


2, 


,043+0-?36 

—0.187 


1, 


■823i}:^i3 


095902.76 


+ 021906.3 


8009 


5.24+0-?5« 

^ — 0.163 


1 


.686+0-049 

^^^—0.048 


0, 


000+0-340 
■'^'-"^-0.000 


095902.76 


+ 021906.3 


8015 


^^^ — 0.183 


1, 


0-1+0.074 
— 0.050 


0, 


091 +0.367 
■^^^-0.021 


095924.46 


+ 015954.3 


8019 


9 O+0.253 
— 0.129 


2, 


.265+0-119 

^^—0.080 


0, 


c: 79 + 1. 144 
■"^'^-0.572 


095924.46 


+ 015954.3 


8020 


2.68+°-}« 

— 0.114 


2, 


.070+0-077 

—0.073 


0, 


000+0-946 
■^""-0.000 


095924.46 


+ 015954.3 


8025 


3.67+0-?09 

— 0.165 


2, 


901+0.088 
■ — 0.082 


0, 


000+0-199 
■'^'-"^-0.000 


095924.46 


+ 015954.3 


8026 


2.98+0-?55 

— 0.144 


2, 


161+0-112 

0.079 


0, 


■88oiJ;gg(, 


105316.75 


+ 573550.8 


1683 


5.42+0-8« 

— 0.713 


2, 


qqi +0.354 
■ — 0.222 


1, 


0.17+3. 163 
■"•^'-1.837 


lUoolO. i 


+ / OOOU.O 


1684 


o r:2 + l-16 
'^■'-^— 0.449 


1 


.587+0-258 

—0.188 


0, 


f;f;o + 1.326 

■""■'-0.663 


105316.75 


+ 573550.8 


4936 




1 


.762+0-024 

^— 0.025 


0, 


000+0-073 
■'^'-"J-0.000 


114656.73 


+ 472755.6 


767 


5.59+°-29? 

— 0.237 


2, 


.068+0-083 

^^^—0.078 


0, 
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Table 4 — Continued 



SDSS Name Obsid Norm T Edge Depth (r) 

(J2000) 10-5 ph keV-i cm'^ g-i at 1 keV 
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Fig. 1. — Histogram showing the distribution of the number of Chandra observations per 
source, for our sample. The axis is logarithmic to clearly illustrate the number of sources 
with more than 2 observations. 
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Fig. 2. — Light curves for sources having 15 or more observations in our sample. 
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Redshift 

Fig. 3. — Redshift and maximum rest-frame time between observations (across all observing 
epochs) for our full sample of Chandra sources. Sources observed in 2 epochs are shown as 
black squares. Those with 3 epochs are plotted in red, and those with 4 or more epochs are 
plotted in green. The solid line at the top represents the rest-frame time from the start of 
Chandra observations (14 Aug 1999 to 1 Jul 2011), which is an upper limit on time scales 
available in the Chandra archive. The second solid line, slightly lower, represents cut-off 
time for an observation to be publicly available for our study. 
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Fig. 4. — 1^2500 as a function of redshift for our full sample of Chandra sources. The symbols 
are the same as in Figure [31 
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Fig. 5. — Black hole ma sses, Mrh^, as a funct ion of redshift. Open circles represent Seyfert 
AGN and quasa r s from iPeterson et al.l (120041 ) , while filled squares repre sent masses taken 
from lShen et al.l (120081 ) for quasars in our sample. Mbh was estimated by lShen et al.l ( 120081 ) 
using the H,^ line for quasars at ^ < 0.7, Mg IIA2800A for 0.7 < z < 1.9, and C IVA1549A 
for z > 1.9. Vertical, dotted lines are plotted at z = 0.7, 1.9 to distinguish the subsamples. 
Some well-known local AGN are identified with star symbols on the plot. 



- 59 - 




Redshift (z) 

Fig. 6. — Median luminosities (estimated as described in §2.4p and redshifts for sources in 
our full sample (black squares) compared to the AOO (red squares), Pao04 (green squares), 
and Pap04 (blue squares) samples. The ?/-axis is logarithmic, and long error bars stretching 
downward off the plot indicate sources that were not detected at high significance. Well- 
known, local AGN (at low redshifts) are plotted as stars using the same colors as in Figure [5l 
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sys 

Fig. 7. — Similar to Figure |6l but showing luminosity as a function of the maximum rest- 
frame time between epochs for each source. We have omitted error bars on luminosity (shown 
in Figure E]) for clarity. For the AOO, Pao04, and Pap08 samples, we determined max(Atsys) 
using observed-fiame times of 14 d, 435 d, and 52 d to span their observation periods; in 
practice, they may have binned individual sources to a smaller time scale resolution. The gray 
shaded region shows the region of time between power spectrum break time scales for AGN 



with {Mbhi bolometric luminosity) = (lO^Mg, 10 



44 



estimated using the relation given in iMcHardy et al.l ( l2006l ) . 



erg s 



and (1O^-^M0, 10^^ erg s 
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2 4 

Redshift z 

Fig. 8. — The mean number of 0.5 — 8 keV counts for sources in Sample HQ as a function 
of redshift. Sources that have been flagged as variable are plotted using red points; non- 
variable sources are plotted as black points. Known BAL quasars are plotted as circles, 
known radio-loud quasars are plotted as triangles, and radio-loud BAL quasars are plotted 
as stars. All other sources are plotted as squares. The solid black hues show typical numbers 
of counts from a hypothetical quasar having an unabsorbed power-law spectrum with F = 2 
in 18 ks of on-axis exposure on an ACIS-I chip, for 0.5-8 keV luminosities Lx = lO'^^, 10^^, 
and 10^^ erg s~^. 
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Mean Counts Per Epoch (HQ) 

Fig. 9. — The solid points show variabihty levels as a function of the number of counts for 
quasars in Sample HQ, using the same symbols as in Figure [HI Sources for which variability 
was detected (not detected) are red (black). The y-axis is the maximum value of \r /tq — 1| 
over all epochs for each source, where r is the observed flux and tq is the best-fit flux. The 
solid black line roughly indicates the 3cr variablity level, which approximates our variability- 
detection criterion. The thick red line indicates the fraction of sources with mean counts 
< X that are identified as variable. 
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Fig. 10. — Each point shows the standard deviation of fractional variation (after accounting 
for the scatter due to measurement errors) in a bin of 100 epoch pairs in Sample HQ as a 
function of the median time between epochs in that bin. The black dot-dashed line shows a 
fit to values at rest-frame times Atgys > 5 x 10''^ s. The red dotted lines indicate la and 3a 
upper limits for 15 pairs of epochs for radio-loud, non-BAL quasars. 
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Fig. 11. — The fractional variation of Sample HQ quasars as a function of time. Red 
and green points indicate known BAL and radio-loud quasars, respectively. We only plot 
points that are > la from zero, in order to clearly demonstrate the amplitudes of fractional 
variation. 
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I F I, Including Error 

Fig. 12. — The x-axis represents the magnitude of fractional variation, F, including measure- 
ment error, observed over time scales Atgys < 5 x 10^ s (black), 5 x 10^ < Atgys < 10^ s (red), 
and Atgys > 10^ s (green). We only consider radio-quiet, non-BAL quasars for this plot. The 
y-axis indicates the upper limit (at 95% confidence) on the true fraction of observations in 
our sample that have \F\ greater than the value on the x-axis. The numbers in the legend in 
parentheses indicate subsample sizes. The upper limits represented by these curves depend 
on several issues, including subsample sizes and measurement errors. 
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Fig. 13. — Excess variance as a function of count-rate luminosity for radio-quiet, non-BAL 
quasars. Sources with excess variance < 0.001 are plotted as open circles y = 0.001. Red 
points are sources at redshift z > 2, while black points are sources at lower redshift. The 
thick green line represents the mean excess variance calculated for bins of 21 sources at 
redshift z < 2. The filled green circle represents the mean excess variance of sources at 
z > 2; the value is less than 0.001. The dashed red lines indicate 2a and 3a upper limits on 
the mean excess variance for quasars aX z > 2. 
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Redshift z 

Fig. 14. — The distribution of Sij = {fij — Cj)/aij for all epochs i of each source j. The 
sources are ordered and binned by redshift. Some points have been clipped off the plot 
boundary for visual clarity. A horizontal dotted (solid) line indicates y = ±1 (y = 0), while 
vertical dotted lines indicate bin boundaries. Red points (connected by red lines) indicate 
the standard deviation of points in each bin, while green points (connected by green lines) 
indicate the rms estimated from the median absolute deviation (MAD). The red points can be 
driven by dramatically variable outliers in a given bin, while the green (MAD) points provide 
a representation of sample variability that is more robust ( §3.5p . Blue lines correspond to 
hypothetical sources that have intrinsic variation of 10%, 20%, or 30% in each epoch. At 
\y\ < 1, measurement error dominates any variability signal in the data. 




52 53 54 

log^Q(Count Luminosity in cts sec-^) 

Fig. 15. — Similar to Figure [TH but showing Sij as a function of count luminosity. 
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Fig. 16. — Change in HR compared to change in luminosity for all epochs of non-BAL, 
radio-quiet quasars in Sample HQ. Squares represent quasars at lower redshifts {z < 2), 
while open circles represent quasars at higher redshifts {z >2). Typical (median) errors are 
shown in the upper right corner. The solid black line represents a linear fit to the full set of 
points, while the red curve represents the two-component toy model described in §3.61 
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Fig. 17. — The change in photon index F between two observations of the same source as 
a function of F in the earher observation. Sources classified as "variable" are plotted in 
red points, while "non-variable" sources are plotted in black. Two outlier points have been 
clipped for visibility. The horizontal dotted line indicates y = and the vertical dotted line 
indicates the median value of F for the data set. All epochs of a single source are plotted 
with the same symbol. 
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Fig. 18. — Same as Figure [HI but showing AF plotted against the change in power law 
normalization, ANorm. 
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Fig. 19. — SDSS spectrum of J123800.91 + 621336.0, plotted in black. The spectrum has 
been smoothed with a boxcar window 11 bins wide. The red spectrum shows the quasar 
composite, normalized to overlap at 5100A. Vertical dotted lines indicate the wavelengths of 
[O III] lines, while Balmer lines are labeled with text. We interpolated over the spectrum at 
4372-4378A due to contamination from sky lines. 
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Fig. 20. — Results of fitting a Galactic- absorbed power law to each epoch of the observed- 
frame (0.5-8 keV) spectra of J1238. Error bars represent la confidence regions. The median 
value of r is shown dotted horizontal line. 
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Fig. 21. — The black histogram shows the distribution of fractional change in counts, 
Fc for segments flagged as (potentially) variable in light curves having mean count rates 
>0.01 counts s~^. The red curve shows the same distribution for our simulated sources 
that have no intrinsic variability. Positive values of Ff. correspond to emission flares, while 
negative values correspond to absorption or dimming events. 
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Fig. 22. — Composite SDSS spectra (plotted in black) for radio-quiet, non-BAL sources with 
redshifts < z < 0.8. The y-axis corresponds to the logarithr a of the composite spectru m, 
normalized to 1 at 5100 A. The SDSS quasar composite from IVanden Berk et al.l (120011 ) is 
overplotted in red for comparison. The top panel shows the composite (in black) for sources 
identified as having at least one variable X-ray epoch. The bottom panel shows the composite 
(in black) for non-variable sources. Vertical dotted lines indicate the wavelengths of [O III] 
lines, while prominent Balmer lines are labeled in the top panel. 
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Fig. 23. — Estimated [O III] equivalent widths as a function of monochromatic luminosity 
at 5100 A for sources classified as variable (red) and non- variable (black). 



